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METHOD AND APPARAITJS F<)R DETERS 

COMPONEPrr USmC CONDmON EWICATOM 

' # * • - 

Background of The Invention 

1. Field of the Invention 

This application relates to the field of vibration analysis and more particulaily to 
performing vibration analysis for the purpose of device monitoring. 

2. Description of Related Art 

The transmission of power to rotors which propel helicopters and other shafts that 
propel devices within the aircraft induce vibrations in the supporting structure. The 
vibrations occur at frequencies that correspond to the shaft rotation rate, medi rat^ bearing 
passing frequency, and harmonics fliereof The vibration is associated with transmission 
OTTor(TE). Increased levels of TE are associated with transmission fail^ Similar types of 
vibrations are produced by transmissdons in fixed installations as welL 

Parts, such as those that may be included in a helicopter transmission, may be 
replaced in accordance! with a predetermined maintenance and parts replacement schedule. 
These schedules provide for replacement of parts prior to failure: The replacement schedules 
may indicate replacement time intervals that are too aggressive resulting in needless 
replacement of working piarts. This may result in incurring unnecessary costs as airplane 
parts are expensive. Additionally, new equipment may have installed faulty or defective 
parts that may feil prematurely. 
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Thus it may be desirable to provide for an efficient technique for detecting part and 
device degradation without unnecessarily replacing parts. It may be desirable that this 
technique also provide for problem detemiination and detection prior to failure. 
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Summary Of The Invention 

In accordance witii one aspect of the invention are a me&od executed in a computer 
system and a computer program product for detemiining a health indicator associated with a 

5 component A plurality of health classifications are determined. At least one condition 
indicator is deternomed quantifying a characteristic of the component A probability 
associated wi& each ofthe health classificatioris is determined The probability is an 
estimation that the component is of a particular health classification given the at least one 
indicator. A determination is fnade as to which of said healdi classifications is associated 

10 with said conq)onent usiog said probabilities associated with said health classifications for a 
given set of observed values. 

In accordance with another aspect of the invention are a method executed in a 
computer system and a computer program product for determining a health status of a 
15 component A plmaHty of condition indicators are selected having a valiie and each having a 
con^onding weightiiig factor, and at least one thrediold val^^ 

classifications. A contribution to a health indicator is determined for each of said condition 
indicators, wherein said determining further comprises, for each of said plurality of 
indicators: detemiining which ofsaid at least two classificatioiis said value of said each 
20 indicator belongs; and determining said contribution to said health indicator by said each 
condition indicator in accordance with a selected one of said at least two classifications and 
said weighted value. The health indicator is determined in accordance witii all contributions 
by each of said condition indicator values. 
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In accordance with one aspect of the invention are a method executed in a computer 
system and a computer program product for detennining an health indicator of a conq>oneiit 
at a subsequent time. A first health indicator of said component at a lime, n, is determined in 
accordance with at least one corresponding condition indicator. A second health indicator of 
tbe component is detennined using a three state Kahnan filter at a time subsequent to time n. 

In accordance with one aspect of the invention are a method executed in a conq>uter 
system and a computer program product for ranking condition indicatora used in detennining 
a health indicator for a component A first set of a plurality of said conditioa indicator is 
detennined. A covariance matrix corresponding to said plurality of condition indicatois is 
detennined. A transfonnation matrix that whitens the covariance matrix is detenniied. 
Differences betwe«i said first plurality of condition indicators and expected values for said 
condition indicators belonging to a health class are detennined using tiie whitening matrix. 
Each health class has a corresponding healtii indicator. A portion of said phirality of 
condition indicators is selected in accordance with those condition indicators haVe the 
smallest of said differences. 

In accordance with one aspect of the invention are a method executed in a computer 
system and computer program product for estimating a conditional indicator value associated 
with a gear pair, the gear pair is modeled as a damped spring model having a contact line 
between said gears. A force, P, is detennined at a point of contact along said contact line 
causing linear and torsional response to each of said two gears in said gear pair. A relative 
movement, d, is determined of said gear pair, in accordance with said force, P, as a sum of 
four responses and a contact deflection, said relative movement d representing a gear model 
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having two degrees of freedom. The relative naovement, d, is used in determining said 
conditional indicator value for transmission error associated witti said gear pair. 

In accordance with another aspect of the invention are a inethod executed in a 
5 computer system and a computer program product for estimating a condition indicator 
associated with a bearing. A bearmg jS^equency ratio is determined for flie bearing. A 
periodic impulse is determined in accordance with the bearing frequency ratio. An intensity 
of an impulse on a bearing surface as a iunctipn of an angle relative to «a bearing fault is 
determined. A decay of a unit impulse is determined, A movement of die bearing is 
determined. A conditional indicator value is determined in accordance with ttie movement 

Li accordance with one aspect of the invention are a method executed in a coniputer 
system and computer program product Tor normalizing a set of at least one observed 
condition indicator. A plurality of conditional indicators and at least one associated factor are 
determined in accordance with previous data acquisitions. A mean and at least one model 
coefScient corresponding to said at least one associated factor are determined. The set of 
obsierved condition indicators is adjusted in accordance with model coefScients and said at 
least one associated factor producing a normalized set of condition indicators. 

In accordance with one aspect of flie invention are a method executed in a computer 
system and computer prograin product for detemiining a condition indicator about a 
characteristic of a component. A distribution of observed data associated with said 
componejit is determined. A difference between said distribution and a normal distribution is 
determined. The condition indicator is determined using the dijSetence. 

.02095633Aaj-> 
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In accordance with another aspect of the invention are a meAod executed in a 
computer system and computer program product for determining a condition indicator 
associated with a component A total impulse signal is determined in accordanpe wifli 
configuration data. The total impulse signal is a supeiposition of gear and bearing noise 
5 represented as a convolution of a gear and bearing signal with a gearbox transfer function. A 
condition indicator is determined in accordance with the total inqtulse sigoaL 

In accordance with yet ano&er aspect of the invention is a niethod executed in a 
computer system and a computer program product for determining a health'Status of a 
10 component using at least one condition indicator. At least one condition indicator is 
determined using at least one of: an impulse determination technique and a statistical 
nonnaUty test The health indicator is detenziined in accordance with the at least ox^ 
indicator. ' 
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Brief Description Of Drawings 

Features and advantages of the present invention will become more qjparent fiom tbe 
following detailed description of exemplary embodiments thereof taken in conjunction with 
the accompanying drawings in which: 

Figure 1 is an exanaple of an embodiment of a system that may be used in performing 
vibration analysis and performing associated monitoring functions; 

Figure 2 is an example representation of a data structure that include aircraft 
mechanical data; 

Figure 3 is an example of parameters that may be included in the type-specific d^ 
portions when the descriptor type is an indexer; 

Figure 4 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is an accelerometer; 

Figure 5 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is a shaft; 

Figure 6 is an example of parameters that may be included in flie type-specific data 
portions when the descriptor type is for a gear; 

Figure 7 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is a planetary type; 
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Figure 8 is an example of parameters that may be included in the type-specific data 
portions when the descriptor type is bearing type; 

Figure 9 is an example of a data structure that includes analysis infomiatiQi^ 

^ ■ ■ * ' 

Figure 1 0 is a more detailed example of an embodiment of a header descriptor of 
Figure 9; 

Figure 11 is an example of a descriptor that may be included in the acquisition 
descriptor groiq> of Figure 9; 

Figure 12 is an example of a descriptor that may be included in the accelerometer 
groiqj of Figure 9; 

Figure 13 is an example of a descriptor that may be included in the shaft desaiptor 
groizp of Figure 9; 

Figure 14 is an example of a descriptor that may be included in the signal average 
descriptor group of Figure 9; ] „ 

Figure 15 is an example of a descriptor that may be included in the envelope 
descriptor group of Figure 9; 

Figure 16 is an example of a planetary gear arrangement; 

8 
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Figure 1 7A is an example of an embodiment of a beaiing; 

Figure 17B is an example of a cut along a line of Figure 17A; 

Figure 1 8A is an example of a representation of data flow in vector transformations; 

Figure 1 8B is an example of a representation of some of the d algorithms that may 
be included in an embodiment, and some of the various ii^uts and ou^uts of each; 

Figure 19 is an exanq>le of a graphical representation of a probability distribution 
function (PDF) of observed data; 

Figure 20 is an example of a graphical representation of a cumulative distribution 
function (CDF) observed data following a gamma (5;20) distribution and the normal CDF; 

Figure 21 is an example of a graphical representation of the difference between the 
two CDFs of Figure 20; 

Figure 22 is an example of a graphical representation of the PDF of observed data 
following a Gamma (5^0) distribution and a PDF of the normal distribution; 

Figure 23 is an example of another graphical representation of the two PDFs fix>m 
Figure 22 shown which quantities as intervals rather than continuous lines; 

9 
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Figure 24A is an example of a graphical representation of the difiFerences between the 
two PDFs of observed data and the normally distibuted PDF; 

Figures 24B-24D are examples of a graphical data displays in coimectioii with a 
healthy S3^em; 

^ ■ ■.' 

Figures 24E-24G are examples of graphical data displays in connection wi& a system 
having a fault; 

Figure 25 is a flowchart of steps of one embodiment for detennining health indicatois 

(His); 

Figure 26 is a graphical illustration of flie probability of a false alarm (PFA) in one 
exanq)le; 

Figure 27 is a graphical illustration of the probability of detection (PD) in one 
example; 

Figure 28 is a graphical illustration of the relationship between PD and PFA and 
threshold values in one embodiment; 

Figure 29 is an graphical illustration of the probability of Ho and threshold values in 
one embodiment; 

Figure 30 is an example of an embodiment of a gear model; 
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Figure 31 is a graphical representation of an estimated signal having an inner bearbig 
fault; and . 

Figure 32 is a graphical representation of the signal of Figure 31 as a firequeiK^y 
spectrum. o 
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Detailed Description of the Preferred Embodimeptfe^ 
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Refemiig now to Figure 1, shoWn is an example of an embodiment of a system 10 
that may be used in performing vibration analysis and monitoring of a macbine such as a 
portion of an aircraft The machine being monitored 12 may be a particular element within an 
aircrafL Sensors 14a througlh 14c are located on the machine to gath^ data from one or more 
componoits of the machine. Data may be collected by the sensors 14a througjb 14c and sent 
to a processor or a VPU16 for data gathering and analysis. The VPU16 analyzes and gafheis 
the data from the Sensors 14a throu^ 14c. 

The VPUl 6 may also use other data in performing analysis. For example, the VPU16 
may use collected data 18. One or more of the Algorithms 20 may be used as input into the 
VPUl 6 in connection with analyzing data such as may be gathered from the Sensors 14a 
through 14c. Additionally, conjBguration data 22 may be used by the VPU16 in connection 
with performing an analysis of the data received for example from the Sensors 14a through 
14c. Generally, configuration data may include parameters and the like that may be stored in 
a configuration data file. Each of these will be described in more detail in paragraphs that 
follow. 

The VPUl 6 may use as input the collected data 1 8, one or more of the algorithnas 20, 
and configuration data 22 to determine one or more condition indicators or CIs. In turn, these 
condition indicators may be used in determining health indicators or His that may be stored 
for example in CI and HI storage 28. CIs describe aspects about a particular component that 
may be useful in making a determination about the state or health of a component as maybe 
reflected in an HI depending on one or more CIs. Generally, as will be described in more 
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detail in paragraphs that follow, CIs and His may be used in connection wifli different 
techniques in determining an indication about monitored components such as Machine 12. 
As described in more detail elsewhere herein, the configuration data may include values for 
parameters that may vary in accordance with the type of the component being mqnitoxed. 

It should be noted that the collected data 1 8 may include data collected over a period 
of time from sensors such as 14a through 14c mounted on Machine 12. A user, such as a Pilot 
26, may use a ^ecial.service processor, such as, the PPU24, connected to the Machine 12 to 
obtain diiT^ent types of data such as the CI and HE values 28. 

As described in connection with Figure 1 , the VPUl 6 may receive inputs from 
Sensors 14a through 14c. These sensors may be different types of data gathering monitoring 
equipment including, for example, high resolution accelerometeirs and index sensors 
(indexors) or tachometers that may be moimted on a component of Machine 12 at carefully 
selected locations throughout an aircraft. Data from these sensors may be sampled at high 
rates, for example, iqp to 100 kilohertz, in order for the VPUl 6 to produce the necessary CI 
and HI indicators. Data from these sensors and accelerometers may be acquired 
synchronously at precise intervals in measuring vibration and rotational speeds. 

Generally, the different types of data gathering equipment such as 14a-14c may be 
sensors or tachometers and accelerometers. Accelerometers may provide instantaneous 
acceleration data along whatever axis on which they are mounted of a particular device. 
Accelerometers may be used in gathering vibration analysis data and accordingly may be 
positioned to optimally monitor vibration generated by one or more mechanical components 
such as gears, shafts, bearings or planetary systems. Each component being monitored may 
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generally be monitored using two independent sensors to provide confinnation of component 
faults and to enable detection of SCTsor faults. 



No accelerometer is completely isolated from any other componemt Thus, the 
component rotational frequencies share as few common divisors as possible in oidff to 
maximize the effectiveness of the monitoring function being perfonned. For example, all 
gears being monitored should have differing number of teeth and aU bearing should have 
differing numbers and sizes of balls or rollers. This may allow individual components to be 
spectrally isolated from each other to the extent that their rotational frequeiides arc'unique. 
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The indexers (mdex sensors) or tachometers may also be used as a particular 
monitoring component 14a througji 14c to gather data about a particular component of 
Machine 12. The indexers produce a periodic analog signal whose frequency is an integer 
multiple of the instantaneous rotation frequency of the shaft that they arc monitoring. These 
15 signals may be generated magnetically using one or more evenly spaced metallic protrusions • 
on the shaft passing by the fixed sensor. Alternatively, these may be monitored optically 
using a piece of optically reflective material affixed to the shaft It should be noted that each 
index point should be fixed in time as precisely as possible. In connection wife magnetic 
sensors, this may be accomplished for example by interpolating the zero crossing times of 
20 each index pulse and similarly for optical sensors by locating either rising or falling edges. 
Assuming the minimal play or strain in the drive train when something is under load, the 
relative position and rate of any component may be calculated using a single index or wave 
form. 
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Because of the high data rates and lengthy processing intervals, diagnostics may be 
performed, for example, on pilot command or oh a predetermined flight regime or time 
intovaL ,^ 

Each of the algorithms 20 produces one or more ds described else^ere hereiii in 
more detail. Generally, the CI may yield useful infomiation about the health of a monitored 
component This condition indicator or CI as well as HI may be used in detemiinin^ or 
predicting faults of dijBferent components. 

It should be noted that the VPUl 6 is intended to be used in a wide variety of 
mechanical and electrical epvironments. As described herein, different component of an 
airoraft may be monitored. However, this is only one example of a type of environment in 
which the system described herein may be used. As known to fliose skilled in ihc art, the 
general prindples and techniques described herein have much broader and general 
applicability beyond a specific aircraft environment that may used in an exanQ)le here. 

In connection wifli the use of CIs, the VPUl 6 uses the ds as input and portions of the 
data such as, for example, used in connection with an algorithm to provide His. These are 
described in more detail in paragr^hs that follow. 

It should be noted that in a particular embodiment, each mechanical part being 
mom'tored may have one or more sensors associated with it where a sensor may include for 
example an accelerometer or a tachometer. Generally, accelerometers may be used, for 
example, to obtain data regarding vibrations and a tachometer may be used, for example, to 
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gain information and data regarding rotation or speed of a particular object. Data may be 
obtained and converted Srom the time to the frequency domain. 

A particular algorithm may provide one or more ds. Each of the algorithms may 
5 produce or be associated with a particular CL One or more (Us may be used in combination 
with a function to produce an HI for a particular part or type. As will be described in more 
detail herein, each of the algorithms may be associated or classified with a particular part or 
type. The d generally measures vibrations and applies a function as described in accordance 
for each algorithm. Generally, vibration is a function of the rotational firequency in the 
10 amount of torque. Using torque and a particular frequency, a CI is appropriately determined 
in accordance with a selected algorithm for a part 

The algorithms 20 may be classified into four families or groups in accordance with 
the different types of parts. In this example, the families of algorithms may include shaft, 

15 gears, bearings, and planetary gears. Associated with each particular part being monitored 
may be a number of ds. Each CI may be the result or output of applying a different one of 
the algorithms for a particular family. For example, in one embodiment, each gear may have 
an associated 27 CIs, each bearing may have 19 CIs, each shaft may have 22 ds, and each 
planetary gear may have two or three ds. It should be noted that each one of these numbers 

20 represents in this example a maximirai number of CIs that may be used or associated with a 
particular type in accordance with the number of algorithms associated wdth a particular class 
or family. Generally, the different number of CIs that may be associated wifli ai particular 
type such as a gear try to take into account the many different ways in which a particular gear 
may fail. Thus, a CI reflects a particular aspect or characteristic about a gear with regard to 
25 how it may fail. 
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Different techniques used in computing CIs are described, for example, in "Introduction to 
Machinery Analysis and Monitoring, Second Edition", 1993, Perm Well Publishing Company 
of Tulsa, OK, ISBN 0-87814-401-3, and "Machinery Vibration: measurement and analysis", 
1991, McGraw-Hill Publishing, ISBN-0-07-07193^5. 

Referring now to Figure 2, shown is an example of a data structure 50 that includes 

aircraft mechanical data. Generally, this data structure includes one or more descriptors 56a 

> . > ■' ' ■■ ... 

through 56n. In this embodiment there may be one descriptor for -each sensor. A descriptor 

associated with a particular sensor includes the parameters relevant to ihc particular 

10 component being monitored. Each of the descriptors such as 56a includes three portions of 

data. The field 52 identifies a particular type of descriptor. Each of the descriptors also 

includes a cornmon data portion 54 which includes ihose diata fields common to all descriptor 

types. Also included is a type specific data portion 56 which includes different data fields, 

for example, that may vary in accordance with the descriptor type 52. . 

15 

Descriptor types may include, for example, an indexer, an accelerometor, a shaft, a 
gear, a planetary gear, or a bearing descriptor type value corresponding to each of the 
different types of descriptors. The conGunon data portion 54 may include, for example, a 
name, part number and identifier. In this example, the identifier in the common data filed 54 
20 may uniquely identify the component and type. 

Referring now to Figures 3 througih 8, what will be described are examples of 
descriptor type specific parameters or infonnation that may be included in a descriptor of a 
particular type, such as in area 56 of the data structure 50. 

25 
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Referring now to Figure 3, shown is an example of parameters that may be included 
in a descriptor 60 which is an indexer descriptor type. The parameters that may be included 
are a channel 62, a type 64, a shaft identifier 66, a pulses per revolution parameter 68, a pulse 
width parameter 70, and a frequency of interest 72 for this particular type of descriptor^ It 

5 should be noted that the type in this example for the index or descriptor may be one of 
sinusoidal, pulse such as 1/ rev, or optical. The shaft identifier 66,is that as may be read or 
viewed by the indexer that calculates the shaft rate. The pulse width 70 is in seconds as the 
unit value. Additionally, the fi-equencyofinterest 72 for this descriptor type is a nominal 
pulse frequency that is used in computing the data quality signal to noise ratio. The use of 
10 these particular data stractures and parameters is described in more detail in paragr^hs that 
follow. 

Referring now to Figure 4, shown is an exanq)le of tiie parametCTS that may be 
included in an accelerometer descriptor type 80. The descriptor for an accelerometcr type 
15 may include the channel 82, a type 84, a sensitivity 86 and a frequency of interest 88. In this 
example for the accelerometer descriptor type, the type may be one of normal, or remote 
charge coupled. The frequency of interest may be used in computing the data quality signal 
to noise ratio. The frequency of interest for a gear is the mesh rate which may be cdctdated 
fi-om the gear shaft rate and the number of teefli of the gear. 

20 



Referring now to Figure 5, shown is an example of descriptor type specific parameters 
or data that may be included when a descriptor type is the shaft descriptor. A shaft descriptor 
90 includes path parameter or data 92 and nominal RPM data 94. The path data is an even 
length sequence of gear tooth counts in the mechanical path between the shaft in question and 
25 a reference shaft. The driving gears alternate with driven gears such that the expected 
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frequency of a gear, shaft, bearing and the like may be determined based on an input shaft 
RPM. 



Referring now to Figure 6, shown is an example of data of parameters that may be 
included in a descriptor when the descriptor type is the gear descriptor. Included in the gear 
descriptor 100 is the shaft identifier 102 to which the gear is mounted and a parameter 104 
indicating the number of teeth in the gear. 

Referring now to Figure 7, shown is an example of an embodiment of a planetary 
descriptor 110 identifying those parameters or data that may be included when the type is a 
planetary descriptor type. The planetary descriptor 1 10 may include an input shaft identifier 
1 12, an output shaft identifier 1 14, a parameter indicating the number of planet gearsLllfi, a 
parameter indicating the number of teeth on the planet gear, a parameter 120 indicating the 
number of teeth on the ring gear, and a parameter 122 indicating the number of teeth on the 
sun gear. It should be noted that the number of teeth on a planet gear relates to a planet 
carrier that is assumed to be mounted to the output shaft. Additionally, the ring gear is 
described by parameter 120 is assumed to be stationery and the sun gear 122 as relied to 
parameter 122 is assimied to be mounted to the input shaft. It should be noted that the path 
between the input and the output shaft may be reduced to using a value S for the driving path 
tooth count and R+S as the driven path tooth coimt where R and S are the ring and sim tooth 
coimts respectively. An example of a planetary type gear is described in more detail 
elsewhere herein. 

Referring now to Figure 8, shown is an example of a bearing descriptor 130. The 
bearing descriptor 130 may include descriptor type specific fields including a shaft identifier 
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132, a cage ratio 134, a ball spin ratio 136, an outer race ratio 138 and an inner race ratio 140. 
An example of a bearing is described in more detail elsewhere herein. 

r 

It should be noted that the data structures described in connection wifli Figures 2 
through 8 are those that may be used in storing data obtained and gathered by a sensor such * * 
as 14a when monitoring a particular component of a machine 12. 'Data maybe gathered and 
stored in the data structure for a particular descriptor or descriptors and sent to the VPU 16 
for processing. It should be noted that a particular set of data may be gathered at a particular 
instance and time, for example, in connection with the synchronous data ga&eiing described 
elsewhere herein. In connection with this, a data set may include multiple descriptois from 
sampling data at a particular point in time which is sent to the VPU 16. 

What will now be described are those data structures that may be associated with an 
analysis definition that consists of a specific data acquisition and a subsequent processing of 
this data to produce a set of indicators for each of the desired components. 

Referring now to Figure 9, shown is an example of the data stmcture 150 that contains 
analysis data. Each instance of analysis data 150 as represented ia the data stracture includes 
a header descriptor 152 and descriptor groups noted as 164. In this example fhere are five 
descriptor groups althougji the particular nmnber may vary in an embodiment Each of the 
descriptor groups 154 through 162 as identified by the group identifier 164 includes one or 
more descriptors associated with a particular group type. For example, descriptor group 154 
is the acquisition group that includes a descriptor for each sensor to be acquired. The 
accelerometer group 1 56 consists of a descriptor for each accelerometer to be processed. The 
shaft group 158 includes a descriptor for each shaft to be processed. The signal average 
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group 1 60 includes a descriptor for each unique parameter set The envelope ffoup 162 
includes a descriptor for each unique paraineter. 

Referring now to Figure 10, shown is a more detailed example of a header descriptor 
170. Parameters that may be included in a header descriptor 170 include: an analysis 
identifier 172, acquisition time out parameter 174 and processing ^time out paranoeter 176. Jn 
this example, the acquisition, time out and processing time out parameters are in seconds. 

Refening now to Figure 1 1 , shown is an example of a descriptor that may be included 
in the acquisition group. A descriptor 180 included in the acquisition groxtp may include a 
sensor identifier 182, a sanq^le rate parameter in Hz 184, a sample duration in seconds 186, a 
gain control setting, such as "auto" or "fixed" 1 88, an automatic gain control (AGQ 
acquisition time in seconds 190, an automatic gain control (AGC) headroom factor as a 
number of bits 192 and a DC offset compensation enable 194. 

Referring now to Figure 12, shown is an example of a descriptor 200 that may be 
included in the accelerometer group. A descriptor in the accelerometer group may include a 
parameter that is an accelerometer acquisition analysis group identifier 202, a list of 
associated planetary identifiers to be processed 204, a list of associated shaft an^ysis groiqi 
identifiers to be processed 206, a processor identifier 208, a transient detection block size 
210, a transient detection RMS factor 212, a power spectrum decimation factor 214 specified 
as a power of 2 and a power spectrum block size also si>ecified as a power of 2. 
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DTohe oSbodiniOTt, the lisTof associated planetary identifier 204 also includes two 
signal average analysis group identifiers for each planetary identifier, first identifier 
corresponding to the input shaft and a second corresponding to an output shaft 

5 It should be noted that the processor identifier 208 will be used in connection with 

assigning processing to a particular DSP or digital signal processor. 

Referring now to Figure 13, shown is an example of an embodiment of a descriptor 
280 that may be included in the shaft group. The descriptor 220 may include a shaft 

) identifier 222, a signal average analysis group identifier 224, a list ofgear identifiers to be 
processed 226, a list of bearing identifiers to be processed 228 and a Ust of associated 
envelope analysis group identifiera 230. 

Referring now to Figure 14, shown is an example of a descriptor 232 that may be 
included in the signal average group. It should be noted that the signal average group 
includes a descriptor for each unique parameter set The signal average processing group is 
run for each accelerometer and shaft combination even if it has the same parametera as 
another combination. Each descriptor 232 may include a number of output points per 
revolution 234 and a number of revolutions to average 236. 

Referring now to Figure 1 5, shown is an example of a descriptor 240 that may 
included in the envelope group. It should be noted that the envelope group includes a 
descriptor for each umque parameter. It is not necessary to repeat an envelope processing for 
each bearing if the parameters are the same. Each descriptor 240 may include a duration 
parameter 242 specifying the seconds of raw data to process, an EFI size 244 which is a 



22 



WO02/695«3 PCTAJS02/16380 
power of 2, a lower bound frequency in Hz 246, and an upper bound finequency, also in, Hz 
248. .... 

Referring now to Figure 16, shown is an example of an embodiment 300 of a 
5 planetary gear arrangement Generally, a planetary gear arrangement as d^cribed in 

connection with the dififerent types of gears and items to be monitored by the system 10 of 
Figure 1 may include a plurality of gears as configured, for example, in flie embodimCTt 300. 
hicluded in the arrangement 300 is a ring gear 2102 a plurality of planet gears 304a through 
304c and of sxm gear 306. Generally, the gears that are designated as planets move around 
10 the sun gear similar to that as a solar system, hence the name of planet gear versus sim gear. 
The arrangement shown in Figure 16 is a downward view representing the different tjpes of 
gears included in an arrangement 300. 

Referring now to Figure 1 7 A, shown is an example of an embodiment 320 of a 
15 bearing. The bearing 320 includes a ring or track having one or more sphqrical or cylindrical 
elements (rolling elements) 324 moving in the direction of circular rotation as indicated by 
the aiTOWs. Different characteristics about such a structure of a bearing maybe important as 
described in connection with this embodiment One characteristic is an "inner race" which 
represents the circvunference of circle 322a of the iimer portion of the ring. Similarly, the 

- 20 "outer-race~or circumference 322b representing the^uter portion of the ring may be a 

consideration in connection with a bearing. 

Referring now to Figure 1 7B, shown is an example of a cut along line 17B of Figure 
17A. Generally, this is cut through the ring or track within which a bearing or bearings 324 
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rotate in a circular direction. The ball bearings move in unison with respect to the shaft 
. within a cage that follows a track as well as rotate around each of their own axis. 

Referring now to Figure 1 8A, shown is an example of a representation 550 of 
5 different transfoimations that may be perfoimed and the associated data flow and 

dependencies for each particular sensor. The output of the transformations are transfonnatioii 
vectors and may be used in addition to analysis data or raw data, such as bearing fi^ency, 
mesh frequency, and the like, by an algorithm in producing a CL 

Referring to the representation 550, an in going arrow represents data flow input to a 
10 transformation. For example, the FF or Fast Fourier transfonn takes as an input data fiom the 
Al signal average data transform. Al has as input the accelerometer data AD. It should be 
noted that other embodiments may produce different vectors and organize data inputs/ou^uts 
and intermediate calculations in a variety of different ways as known to those skilled in the 
art. 

15 

Referring now to Figure 18B, shown is an example of a representation 350 relating 
algorithms, a portion of input data, such as some transfomiation vectors, and ds produced 
for each type of component, that may be included in an embodiment. Other embodiments 
may use different data entities in addition to those shown in connection with Figure 18B. As 

20 described elsewhere herein, each type ofcomponent in this example is one of: indexer. 
accelerometer, shaft, gear, planetary, or bearing. Certain algorithms may be used in 
connection with determining one or more CIs for more than one component type. It should 
be noted that a variety of different algorithms may be used and are known by one of ordinary 
skill in the art, as described elsewhere herein in more detail. The following are examples of 

25 some of the different techmques that may be used in producing CIs. AdditionaUy, Figure 
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1 SB illustrates an example of relationships between some algorithms, a portion of their 
respective inputs and outputs, as well as how the algorithms may be associated with different 
component types. However, it should be noted that this illustration is not all iijclusive of all 

algoritlmcis, all respective inputs and outputs, aiid all component tjpes. 

■ • • • .• . ' '•*.•.■■* 

What will now be described are algorithms and the one orteore CIs produced fliat 
may be included in an embodiment It should be noted that the nmnber and type of ; 
algorithms included may vary in accordance with an embodiment Additionally, it should be 
noted that Figure 1 8B may not include each and every input and output for im algorithm as 
described herein and other embodiments of the algorithms described generally herein may 
also vary. 

The data quality (DQ) algorithm 356 may be used as a quality assurance tool for the 
DTD d. DQ performs an assessment of the raw uncalibrated sensor data to insure that the 
entire system is perforaiing nominally. DQ may be used to identiJ^, for exampl^ bad wiring^ 
connections, faulty sensors, clipping, and other typical data acquisition problems. The DQ 
indicator checks the output of an accelerometer for 'T^ad data". Such "bad data" causes the SI 
to be also be "bad" and should not be used in determining health t:alcidations. 

What will now be described are the different indicators that may be included in an 
embodiment of the DQ algorithm. ADC Bit Use measures the number of ADC bits used in 
the current acquisition. The ADC board is typically a 1 6 bit processor. The log base 2 value 
of the maximxim raw data bit acquired is rounded up to the next highest integer. Channels 
with inadequate dynamic range typically use less than 6 bits to represent the entire dynamic 
range. ADC Sensor Range is the maximum range of the raw acquired data. This range 
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caimot-exceed- fte operational range of the ADC board, and the threshold value of 32500 is 
just below the maximum permissible value of +32767 or -32768 when the absolute value is 
taken. iJynflwuc Range is similar to the ADC Sensor Range, except the indicator reports 
dynamic channel range as a percent rather than a fixed bit number. CTfRp/n^ indicates the 
5 number of observations of clipping in the raw data. For a specific gain value, the raw ADC 
bit values cannot exceed a specific calculated value. Low Frequency Slope (LowFreqSlape) 
audXow Frequency Intercept (lowFreqInt) use the first 10 points of the power spectral 
density calculated fiiom the raw data and perform a simple linear regression to obtain the 
intercept and slope in the fi^quency-amplitude domain. SNR is the signal to noise ratio 
10 observed in each specific data channel. A power spectral density is calculated fi^om the raw 
uncaUbrated vibration data. For each data channel, there are known fi^uencies associated 
with certain components. Examples include, but are not limited to, gear mesh frequences, 
shaft rotation rates, and indexer pulse rates. SNR measures the rise of a known tone 
(corrected for operational speed differences) above the typical minimum baseline levels in a 
15 user-defined bandwidth (generally +/- 8 bins). 

The Statistics (ST) algorithm 360 is associated with producing a plurality of statistical 
indicators 360a. The Root-Mean-Square (RMS) value of the raw vibration amplitu<te 
represents the overall energy level of the vibration. The RMS value can be used to detect 
20 major overall changes in the vibration level. The Peak-To-Peak value of the raw vibrating 
amplitude represents the difference between the two vibration extrema. When failures occur, 
the vibration amplitude tends to increase in both upward and downward directions and thus 
the Peak-To-Peak value increases. The Skewness coefficient (which is the third statistical 
moment) measures the asymmetry of the probability density function (p.d.f.) of the raw 
vibration amphtude. Since it is generally believed that the p.d.f is near Gaussian and has a 
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Skewness coefficient of zero, any large deviations of this value from zero may be an 
indication of faults. A localized defect in a machine usually results in impulsive peaks in the 
raw vibration signal, which affects the tails of the p.d*f. of the vibration amplitude. The fouitii 
moment (Kurtosis) of the distribution has the ability to enhance the sensitivity of such tail 
changes. It has a value of 3 (Gaussian distribution) when flie machinery is healthy. Kurtosis 
values larger than 3.5 are usually an indication of localized defects. However, distributed 
defects such as wear tend to smooth the distribution and thus decrease the Kurtosis values. 

* ' ' .* 

The ST algorithm may be performed on the following vectors: AD raw acceleroineter 
data, Al signal average data, RS residual data, NB narrow band data, and EV envelope data 
and others, some of which are listed in 360b. 

The Tone andBase Energy algorithm(TB) 362 uses tone energy and base energy. 
Tone Energy is calculated -as the sum of all the strong tones in the raw vibration spectrum. 
Localized defects tend to increase the energy levels of the strong tones. This indicator is 
designed to provide an overall indication of localized defects. "Strong tones" are determined 
by applying a threshold which is set based on the mean of all the energy contents in the 
spectrum. Any tones that are above this threshold are attributed to this indicator. The Base 
Energy measures the remaining energy level when all the strong tones are rOTioved frcmi the 
raw vibration spectrum. Certain failures such as wear, do not seem to ^ect the strong tones 
created by shaft rotation and gear mesh, the energy in the base of the si>ectrum could 
potentially be a powerful detection indicator for wear-related failures. Note that the sum of 
Tone Energy and Base Energy equals the overall energy level in the spectrum. 
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SI are miscellaneous shaft indicators. SOI (Shaft Order 1 in g) is the once-per-rev 
energy in the signal average, and is used to detect shaft imbalance. S02 (Shaft Order 2 in g) 
is the twice-per-rev energy in the signal average, and is used to detect shaft misalignment 
GDF (Gear detector fault) may be an effective detector for distributed gear faults such as 
wear and multiple tooth cracks, and is a complement of the indicator signalAverageLl (also'- 
known as gearLocalFault). » 

In addition to the specifically referenced vectors below, the SI algorithm takes input 
from the indexerzCTo-crossing vector (ZC). • . 

^demodulation analysis (DA4) 370 is designed to farther Tcvea^^ 
modulation by using the Hilbert transform on either the narrow band signal (narrow band 
demodulation) or the'signal average itself (wide band demodulation) to produce the 
Amplitude Modulation (AM) and Phase Modulation (FM) signals. The procedures involved 
to obtain such signals are: 

Perform Hilbert transform on the narrow band signal (or signal average). 

Compute the amplitude of the obtained complex analytic signal to obtain the 
AM signal. 

Compute the phase angles of the analytic signal to obtain the FM signaL 

Compute the instantaneous amplitude of the analytic signal to obtain the dAM 

signal. 

Compute the instantaneous phase angles of the analytic signal to obtain the 
dFM signal. 
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The DM algorithm is performed on the band passed jQltered data at a frequency of interest by 
taking a Hilbert Window function of the frequency domain data and conv^ting the data bade 
to the time domain. • 

The Sideband Modulation (SM) 368 analysis is designed to reveal any sideband 
activities that may be the results of certain gear faults such as eccentricity, misalignment; or 
looseness. * 
CIs included in 368a dicDSMn. DSMn is an indicator that characterizes the Degree of 
Sideband Modulation for the nth sideband (n = 1, 2, and 3). The DSMn is calculated as flie 
sum of both the nth higji and low sideband energies aroimd the strongest gear meshing 
hamionic. As indicated in 368b, the SM algorithm is performed on the Fast Fourier transfomi 
vector (FF). 

The Planetary Analysis (PL) 364 extracts .the Amplitude Modulation (A^ 
produced by individual planet gears and compares the "uniformity" of all the modulation 
signals. 

In general, when each planet gear orbits between the sxm and the ring gears, its vibration 
modulates the vibration generated by the two gears. It is believed that when one of the planet 
gears is faulty, the amplitude modulation of that planet gear would behave differently than 
the rest of the planet gears. The procedure to perform this algorithm is to obtain signal 
averages for the input, output, and planet shafts. For each signal aiverage: 
Locate the strongest gear meshing harmonic. 

Bandpass filter the signal average around this frequency, with the bandwidth 

equals to twice the number of planet gears. 

Hilbert transform the bandpass filtered signal to obtain the AM signal. 
Find the maximum(MAX) and minimum(MIN) of the AM signal. 
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Calculate the P7fl77er Ge^ir FflwA (PGi9 indicator as inclu^^ 
according to the equation PGF = MAX(AM) / MIN (AM). 
The inputs to the PL algorithm are the raw accelerometer data (AD) and the indexer zero- 
crossing data (ZC). 

The Zero-Crossing Indicators (ZI) algorithm 354 is performed on the zero-crosising ' 
vector (ZC). The zero crossing indicators may be determined as follows: 
Dj^In y^, ^Inj^j=^ 0. JV — 2 , the stored zero-crossing intervals 

pulseIntervalMean^Mean{D) , • 

The Shaft Indicators (SI) algorithm 358 calculates miscellaneous shaft indicators 
included in 358a. SOI (Shaft Order 1 in g) is the once-per-rev energy in the signal average, 
and is used to detect shaft imbalance. S02 (Shaft Order 2 in g) is the twice-per-rev energy in 
the signal average, and is used to detect shaft misalignment 

SOS (Shaft Order 3J, is the three-per-rev energy in the signal average, and is used to . 
detect shaft misalignment. The miscellaneous shaft indicators may also be included in an 
embodiment defined as follows: 

p ^ nuinPathPairs 

Y\shaftPath^f 

shaftRatio^^ ^^^^^S 

— ^ f~i - driven 

TlindexPath2i 
mdexRatio = -ff — 

, . jy , indexRatio , ^ 

dnveRotioss . pulsesUsed 

shaftJiatio 



30 



WP 02/095633 PCTAJS02/16380 

60 



shaftSpeed = 



pulselntervalMean * driveRatio 



resampleRaie = ^^^fi^P^^ , pointsPerRev 
60 

i?^= residual data, 

5 Al = signal average, 

signalAverageLl = ^P{^} 
Rms{Al) 

/y= FFT of the signal average, 
shqftOrderj^^FFj , 

iO gearDistFault^ ^^^^f^} 

Stdev{AL) 



As described elsewhere herein, ^eari)iy£Ffl«ft (t?^^ 
distributed gear faults such as wear and multiple tooth cracks, and is a complement of the 
15 indicator signaL^verageLl (also knovm as gearIx)caIFaidt), 

In addition to the specifically referenced vectors below, the SI algorithm takes iiq>ut 
from the index er zero-crossing vector (ZC) and may also use others and indicated above. 

20 The following definitions for indicators may also be included in an embodiment in 

connection with the SI algorithm: 

shaftPath is defined for the shaft descriptor 

indexPaih is the path of the shaft seen by the indexer used for signal averaging 
numPaihPairs is the number of path pairs defined for shaftPath ^dindexPath 
25 /?w&e5'L&e£/ is the number of pulses used pCT revolution of the indexer shaft 

pulselntervalMean is the mean of the zero-crossing (ZC) intervals 
pointsPerRev is the number of output points per revolution in the signal average. 
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The Bearijig Energy (BE) algorithm 376 performs an analysis to reveal the four 
beaiing defect frequencies (cage, ball spin, outer race, and inner race frequencies) that usually 
modulate the bearing shaft frequency. As such, these four frequencies are calculated based on 
the measured shaft speed and bearing geometry. Alternatively, theibur frequency ratios may 
be obtained from the bearing manufacturers. The energy levels associated with these four 
frequencies and their harinonics are calculated for bearing fault detection. They are: 

- Cage Energy: the total energy associated with the bearing cage defect 
frequency and its harmonics. Usually it is detectable only at the later stage 
of a bearing failure^but some studies show that this indicator may increase 
before the others. 

- Ball Energy: the total energy associated with the bearing ball spin defect 
frequency and its harmonics. 

- Ow/erifcce-E/jergy.rthe total energy associated with the bearing outer 
defect frequCTcy and its harmonics. 

- Inner Race Energy, the total energy associated with the bearing inn^ race 
defect-frequency and its harmonics^- 

The Total Energy indicator gives an overall measure of the bearing defect energies. 

In one embodiment, one or more algorithms may be used in determining a CI 
representing a score quantifying a difference between observed or actual test distribution data 
and a normal probability distribution function (PDF) or a normal cumulative distribution 
function (CDF). These one or more algorithms may be categorized as belonging to a class of 
algorithms producing CIs using hypothesis tests ("hypothesis testing algorithms") that 
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provide a measmre of difTerence in determining whether a given distribution is not normally 
distributed. These h^'pothesis testing algorithms produce a score that is used as a CI. The 
score may be described as a sum of differences between an observed or actual test 
distribution function based on observed data and a normal PDF or normal CDF. An 
5 algorithm may exist, for example, based on each of the following tests: Chi-Squared 
Goodness of fit (CS), Kolmogorov-Smimov Goodness of fit (KS), lilliefors test of 
normality, and Jarque-Bera test of normality (JB). Other embodiments may also include 
other algorithms base<^ on other tests for normality, as known to those of ordinary skill in the 
art The hypothesis tests compare the test distribution to the normal PDF, for example as 
10 with CS test, or the normal CDF, for example as with the KS and Lilliefor tests. 

■ ^ ■ 

What wall now be described is an example in which the CS test is used in determining 
a score with a test distribution of observed actual data. In this example, the test distribution 
of observed data forms a Gamma (5, 20) distribution function, having and alpha value of 5 
15 and a beta value of 20. The mean of this Ganmia(5,20) distribution is alpha ♦beta having a 
variance of alpha * beta^ . The Gamma (5,20) distribution fimction is a tailed distribution 
which graphically is similar to that of a normal distribution. 

Refeiring now to Figure 19, shovm is an example of a graphical r^resentation 400 of 
20 observed data. 

Referring now to Figure 20, shown is an example of a graphical representation 410 of 
the normal CDF and the Gamma (5,20) CDF of random data. Referring now Figure 21, 
shown is an example of a graphical representation 420 of the difference between the normal 
25 CDF and the Gamma (5,20) CDF. 
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In one embodiment, if there are 1 000 test samples used in forming a single 
CDF, the graphical representation, for example, in Figure 21 represents differences in 1 000 
instances where the difference between the expected value (Normal CDF) and the maximum 
deviation of the (in this case defined as the score) observed gamma CDF can exceed some ' 
critical value. The critical value is that statistic which represents some predefined alpha eiror 
(the probability that the test indicates the distribution is not normal when in fad it is noimal - 
this is typically set at 5%.) If the score exceeds the critical value, the distribution is said to be 
not normal statistic. The score is the maximum deviation fi-om this statistic or alpha' value. 

It should be noted that the sensitivity or goodness of the test increases as the number 
of samples or instances (degrees of fi-eedom "n") increases approximately as the square root 
of "n". For example; in the case where 1000 instances or samples are used such that n=1000, 
the sensitivity or ability of this CI to be used in detecting gear faults, for example, is roughly 
3 1 times more powerful than kurtosis in identi^g a non normal distiibutioiL 

As another example, in the algorithm usmg the CS test, flie normal PDF is used. 
Refemng now to Figure 22, shown is a graphical representation 430 of the normal PDF and 
the PDF of the Gamma (5,20) distribution. The representations of Figure 22 are diawh as 
continuous lines rather than discrete intervals. _ 1. 

Referring now to Figure 23, the quanti'ties of the x-axis represented in Figure 22 are 
shown in another representation 440 as being divided into discrete bins, intervals, or 
categories. For example, there may be 4 bins or intervals between any two integer quantities. 
Between 0 and 1, bin 1 includes values between [0,0.25), bin 2 includes values between 
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[0.25, 0.50), bin 3 includes values between [.050,0.75) and bin 4 includes values between 
[0.75, 1 .0). For each bin,, determine the number of observed and expected values, and their 
difference. Square each of the differences for each bin and then add all the differences and 
divide by the expected value for each bin. The CS test which sums all the differences for eadh 
5 category divided by the expected value for each category represented as: V* 

for k categories or bins, k- 1 degrees of freedom, fi is observed data and ei is expected data 
value or niunber in accordance with a normal distribution. 

For each bin, take the difference between the obsearved and expected observation^ 
Square this value and divided by expected number of observation. Sum over all bins. The 
statistic, the critical ^value is the at k-1 degrees of freedom may be^ for example, 90.72 
which is much greater than the .05 alpha value of a x^ which is 54.57 for 39 degrees of 
freedonl or 40 categories/bins. Thus, the observed data in this example as indicated by the 
statistic is not normally distributed. Figure 24A represents graphically a di£ference between 
observed and expected values for each bin or interval of Figure 23. 

It should be noted that the foregoing algorithms provide a way of measuring both the 
skewness and kurtosis simultaneously by comparing the PDF or of the test distribution 
against-tiie PDF/CDF of a standard normal distribution in which a score is used as a CI as 
described above. 

As known to those of ordinary skill in the art, other algorithms belonging to the 
hypothesis testing class may be used in computing CIs. The particular examples, algorithms, 
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and tests selected for discussion herein are representations of those that may be included in 
the general class. 

What will now be described is another algorithm that may be used in determining a CI 
5 in an embodiment of the system of Figure 1. This may be refoxed to as an impulse 

determination algorithm that produces a CI indicating an amount of vibration that may be 
used in detecting a type of fault The impulse detennination algorithm takes into account the 
physical model of the system. One type of fault that this technique may be used to detect is a 
pit or spall on eithen gear tooth, inner bearing race, outer bearing race or bearing roller 
10 element. This technique uses a model designed to detect this type of fault where &e modd is 
based on knowledge of the physical system. For example, if there is al pit or spall on a 
bearing, this may produce a vibration on a first bearing which may further add vibrations to 
other components connected to or coupled to the bearing. 

IS 

In one embodiment, a model can be detomined for a particular configuration by using 
configuration data, for example. In one configuration, for example, a signal received at a 
sensor may be a superposition of gear and bearing noise fljat may be represented as a 
convolution of gear/bearing noise and a convolution of the Gear/B earing signal witii the 
20 geaibox transfer function. Giverithis, ifonetype offaultisapitorspalloneithera: gear 
tooth, inner bearing race, outer bearing race or bearing roller element, a model that is 
designed to look for this type of fault can take advantage of knowledge of the physical 
system. 

The impulse detennination algorithm uses Linear Predictive Coding (LPQ 
techniques. As known to those skilled in the art, LPC may be characterized as an adaptive 
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ti^je ofsignal processing algorithm used to deconvolute a signd into its base In 
the case of a pit/spall fault, the base signal components are an impulse train generated by the 
fault hitting a surface (e.g gear tooth with geartooth, inner race with roller element, etc) and 
the bearing/case transfer function. The bearing, geiar and case have th^ own transfer 
5 functions. Convolution here is transitive and multiplicative. As such, LPC techniques may 
be used to estimate the total convolution function of the total vibration that may be induced. 

For example, in this arrangement, the total amount of vibration representing the total 
impulse signal generated by a configuration may be represented as: 

[z/npwZse]® f {Gear}® /(Bearing)® f (Case) s [impulse]® {/(Gear)® f (Bearing)® f(Case)] 

in which ® represents the convolution operation. 

It should also be noted that convolution is a homomorphic systetn such that it is 
monotonically increasing and that logarithmic transformations hold. Thus the relationship of 
c = a*b also holds for Log c = Log a + Log b. A "dual nature" of convolution is used in 
following representations to equate operations using convolution in the time domain to 
equivalent multiplication operation in the frequency domain. 

If "j^ represents the to^^^^ tlie 
response of the system for a series of elementary input impulses "imp" such that y is the 
convolution of imp and h, then this may be represented as: 

y = imp®h 
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and then converting "y" and "h" each, respectively, to the frequency domain 
represented as "Y" and "H", as may be represented by the following: 

y = 3(y),jfir = 3(A) 

taking the Fourier transform (FFT) of each where H represents the transfer function. * 
The convolution in the time domain may be equated to a multiplication in the frequency 
domain represented as: 

in which IMP is the Fourier transformation of imp into the frequency domain. Above, wop is 
in the time domain. 

The convolution in the time domain is equivalent to multiplication in the Frequency 
Domain. Referring to the homomorphic property of convolution, it follows /Quit: 

log(y) = Iog(/MP)+ log(/fX 
therefore 

log(iMP) - log(y)- log(ifX 
/il4P = exp(log(r)- ]og(/f )) 

and finally 

Using the foregoing, the system transfer function "H" may be estimated for the 
Gear/Bearing and Case to recover the impulse response allocated vsdth a Gear or Bearing 
pit/spall fault. The estimation of this transfer function "H" maybe accomplished using 
Linear Predictive Coding (LPC) techniques. LPC assumes that the Transfer Function is a 
FIR filter, and as such, the auto-correlation of the time domain signal maybe used to solve 
for the filter coefficients in a minimum sum of square error sense. 
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Using the UPC model, there is an impulse that is convoluted with a FIR filter, such 



that: 



y[n]-a^x[n-l]'^a2x[n-2]^a^x[n-3]'h.., 

LPC techniques may be used to estimate the coefEcierits a=(aj an) for an order pin a 
minimum sum of square error sense, n=p+l. The standard least squares enror estimators may 
Id be used, wherein y = y[l, 2, • . . n], and x is the time delayed signal, in which: 



15 



20 



jc[n - 2, /I — 3, n - p - 1] 



where a = (x^ x)'' x^ y. These values for al an may be Used Avith the following equation: 

•* " ■ . * " - 

y^^ = ax, b — (y-^y had^ and the estimator of enor B is: b- . 

Y may also be expressed as: 
; Y=PFT(y(l,2,..n]) 
in which y[l ..n] are values in the time domain expressed in the frequency domain sts a 
Fourier transform of the time domain values. Y represents current time vector measiirements 
in the frequency domain. 

In terms of a and B, the transfer function H may be estimated and r^resented as a/®, 
(freq. Domain). Note that "a" is a vector of the values al . . . an obtained above. 



The homomoiphic property of convolution as described above may be used to 
25 estimate the impulse as represented in: 
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IMP = expOog(Y).log(H)) IMP Equation 

If there is no fault, the impulse, for example, may be characterized 
as "wliite noise". As the fault progresses, the impulse or the value of H becomes larger. The 
CI is the power spectral density at a bearing passing frequency for a bearing fault, or a mesh 
frequency for a gear fault Other CIs based on the foregoing value may be a "score" of the . 
Lilifers test for normality, or other such test 

In the foregoinjg, a pit or spall may cause a vibration or tapping. Subsequently, other 
elements in contact with the ball bearing may also vibrate exhibiting behavior from this initial 
vibration. Thus, the initial vibration of the pit or spall may cause an impulse spectrum to be 
exhibited by such a component having unusual noise or vibration. 

The value of IMP as may be determined using the IMP Equation above represents the 
impulse function that may be used as a "raw" value and at a given frequency and used as an 
input into an HI determination technique. For example, the IMP at a particular frequency, 
since this the spectrum, determined above may be compared to expected values, such as may 
be obtained from the stored historic data and configuration data. An embodiment may also 
take the power spectrum of this raw impulse spectrum prior to being used, for example, as 
input to an HI calculation where the power spectrum is observed at frequencies of interest, 
such as the inner race frequency. For example, if the impulse frmction is within some 
predetermined threshold amoimt, it may be concluded that there is no fault 

What is shown in the Figure 24B and Figure 24C are relative to a healthy system, 
such as a main gearbox, for example, such as in connection with a planetary race fault of an 
SH-60B U.S. Navy Helicopter built by Silorsky, 
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The Figure 24B representation 700 shows an impulse train, in the frequency domain of 

ft 

the healthy system. 

It should be noted that an embodiment may estimate the transfer functioxi H using 
LPC using different techniques. An embodiment may estimate the transfer function H using 
an autocorrelation technique(AutoLPC). An embodiment may also estimate the Pansier 
function H using a cdvariance technique (CovLt*Q. Use of autocorrelation may use less 
mathematical operations, but require more data than using the covari^nce. Alternatively, use 
of the covariance technique may use more mathematical operations but require less data. As 
the amoimt of available data increases, the autocorrelation LPC result converges to the 
covariance LPC result In one example, data samples are at 1 OOKHz with 64,000 data points 
used with the autocorrelation technique due to the relatively large niunber of data points. 

Figure 24C representation 710 shows the data of 700 from Figure 24B in the time 
domain rather than the frequency domain. 

Figure 24D representation 720 shows the power spectral density of the above figures 
as deconvolved time data of frequency v. dB values in a healthy system. 

The foregoing Figures 24B-24D represent data in a graphical display in coimection 
with a healthy system. Following are three additional graphical displays shown in Figures 
24E-24G in connection with an unhealthy system, such as a starboard ring channel which 
exhibit data that may be expected in connection with a pit or spall fault. 
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Figure 24E, representation 730, illustrates an impulse train as maybe associated with 
an unhealthy system in the time domain. Figure 24F, representation 740, iUustpates a 
graphical display of the impulse train in the frequency domain. ' - 

In Figure 24G, shown is an illustration 740 is a graphical Vepres^tation of the power 
spectrum of the impulse train represented in connection with the other two figures for the 
unhealthy system identified by a period impulse train associated with .an inner race bearing • 
fault. In this example, a spike may be viewed in the graphical display as well as the 
harmonics thereof. 

It should be noted that other algorithms and CIs in addition to those described herein 
may be used in producing CIs used in techniques in connection with His elsewhere herein. 

What will now be described is one embodiment in which these CIs may be used. 
Referring now to Figure 25, shown is a flow chart of steps of one embodiment for 
determining the health of a part as indicated by an HI. At step 502, raw data acquisition is 
performed. This may be, for example, issuing appropriate commands causing the VPU to 
perform a data acquisition. At step 504, the raw data may be adjusted, for example;, in 
accordance with particular configuration information producing analysis data as output. It is 
at step 504, for example, that an embodiment may make adjustments to a raw data item 
acquired as may be related to the particular arrangement of components. At step 506, data 
transformations may be performed using the analysis data and other data, such as raw data. 
The output of the data transformations includes transformation output vectMS. At step 508, 
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CIs are computed using the analysis data and transfonnation vector data as may be specified 
in accordance with each algorithm. At step 510, one or more CIs maybe selected. Particular 
techniques that may be included in an embodiment for selecting particular CIs ^ described 
elsewhere herein in more detail. At step 5 1 2, CIs may be normalized. This st^ is described 
5 m more detail elsewhere herein. At step 514, the selected and npmialized CIs are used in 
determining His. Particular techniques for determining His are d^cribed in more detail 
elsewhere herein. 

Ih an embodiment, due to the lengthy processing times, for example in executing the 
different algorithms described herein, HI computations may not be executed in real time. 
Rather, they may be perfonned, for example, when a command or request is issued, such as 
fix>m a pilot or at predetermined time intervals. 

The hardware and/or software included in each embodiment may vary, in one 
embodiment, data acquisition and/or computations may be performed by one or more digital 
signal processors (DSPs) running at a particular clock speed, such as 40MHz, having a 
predetermined numerical precision, such as 32 bits. The processors may have access to 
shared memory. In one embodiment, sensors maybe multiplexed and data maybe acquired 
in groups, such as 8. Other embodiments may vary the number in^ach group for data 
sampling. The sampling rates and durations within an acquisition group may also vary in an 
embodiment. Data may be placed in the memory accessed by the DSPs on acquisition. In one 
embodiment, the software may be a combination of ADA95 and machine codel Prooessois 
may include the VPU as described herein as well as a DSP chip. 
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What will now be described are techniques for nonnalizing CIs in connection with, 
determining His providing more detailed processing of step 5 1 2 as described in connection 
with flowchart 500. 

Transmission erroT (T.E.) depends upon torque. Additionally, vibration depends upoii 
the frequency response of a gear. As such, the CI, which also depends upon T.E. and 
vibration, is a function (generally linear) of torque and rotor speed (which is frequency), and 
airspeed as this may change the shape of the airframe. Thus, technique? that may be used in 
connection with determining the "health state" or HI of a component may noraialize C3s to 
account for the foregoing since His are detemiined using CIs. 

For each bearing, shaft and gear within a power train, a number of CIs may be 
determined. An embodiment may compare CI values to threshold values, apply a weighting 
factor, and sum the weighted CIs to determine an HI value for a component at a particular 
time. 

Because data acquisitions may be made at different torque (e.g. power setting values, the 
threshold values may be different for each torque value. For example, an embodiment may 
use 4 torque bands, requiring 4 threshold values and weTghts for each CL Additionally, the 
coarseness of the torque bands will result in increased, uncontrolled system variance. 
Alternatively, rather than use multiple threshold values and have an uncontrolled variance, an 
embodiment may use a normalization technique which normalizes the CI for torque and rotor 
RPM (Nr), and airspeed, expressed as a percentage, for example, in which a percentage of 
1 00% is perfect Use of these normalized CIs allows for a reduction of configuration such 
that, for example, only one threshold is used and variance may also be reduced. 
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The normalization technique that will now be described in more detail may be used in 
connection with methods of HI generation, such as the non-linear mapping method and flic ; 
hypothesis testing method of HI generation that are also described in more detajl elsewhere 
herein. 

It should be noted that a deflection in a spring is linearly related to the force applied to 
the spring. The transmission may be similar in certain aspects to a large, complex spring. 
The displacement of a pinion and its corresponding Transmission Error (T£.) is proportional 
to the torque applied. TJE. is a what causes vibration, while the intensity of the vibration is a 
function of the frequency response (Nr), where frequency is a function of RPM. Thu^ 
vibration and the corresponding CI calculated using a data acquisition are approximately 
linearly proportional to torque, Nr„ (over the operating range of interest) and/or airspeed 
although at times there may be a linear torque*Nr interaction effect For exan^le, gear box 
manufacturers miay design a gearbox to have minimum TJE. under load, and a graphical 
representation of TJE. vs. Torque is linear, or at least piece wise linear. It should be noted ^ 
that test data , for example used in connection with a Bell helicopter H-1 loss of lube test, 
shows a relationship between CI and torque suggesting linearity. Additionally, tests show that 
airspeed is also relevant factor. Other embodiments may take into account any one or more 
of these factors as well as apply the techniques described herein to other factors that maybe 
relevant in a particular embodiment or other application although in this example, the f^^tcm 
of torque, airspeed and Nr are taken into account 
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An equation representing a model minimizing the sum of square error of a measured 
CI for a given torque value in a healthy gear box is: 

C/ = ^„ + * T orgue + B^Nr + B, Airspeed + TJE... (Equation 1) 

The order of the model may be determined by statistical significance of the 
coefBcients of Equation 1. Li the previous equation, the TJB. of a "healtjiy" component may 
have, for example, a mean of zero (0) with some expected variance. It should be noted that if 
the model fits well for the lower order. Higher order coefBcients are not required and may 
actuaUy induce error.in some instances. The following example is built ' as a first order 
model, higher orders may be solved by extension of that explained in the first order model. 
This model, written in matrix format is: y = B x where 



y = 



cr 



B = [^o Bu.Bn ] andx = 



1 /i Airspeed! 
1 Nk... Airspeed... 

1 Nm Airspeedn 



Each of the CIs included' in the vector y is a particular recorded value for a CI fiom previous 
data acquisitions, for example, as maybe stored and retrieved fiom the collected data 18. 
Also stored with each occuirence of a a for a data acquisition in an embodiment maybe a 
corresponding value for torque (t), Nr, and Airspeed. These values may also be stored in the 
collected data 18. 

The model coefBcients for B may be estimated by minimizing the sum of square error 
between the measured CI and the model or estimated CI using the observed perforaiance 
data. Solving the foregoing for the unbiased estimator of B = (x^'x)"' . The variance of B 
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is: PhrCB) = ^(b -B)(b -B)*" = cr' (x^x)"' where b is an unbiased estimator of Bi The 

unbiased estin^ator of is s^: s' ^^^^ (y-y)'(y-^) = y'y-*>'''y 

n~p—l n-p-l 

In the vector B from y=xB, coefficient Bo represents the noean of the data set for a 
5 particular component which, for example, may be represented as an ofl&et value. Each of the 
oth^ values Bl . Bn are coefficients multiplied by the coiresponding factors, such as 
airspeed, torque, and Nn 

The foregoing B values or coefficients may be determined at a time other than in real- 
10 time, for example, when flying a plane, and then subsequently stored, along with 

corresponding X infomaation, for example, in the collected data store 18. These stored values 
may be used in determining a normalized CI value for a particular observed instance of a 
Clobs in determining an HL The normalized CI may be represented as: 

Clnonnalaed ^ TJE. = ClobS - (B* X) 

15 where Clobs represents an instanqe of a CI being normalized using previously determined 
and stored B and x values. Threshold values, as may be used, for example, in HI 
determination, may be expressed in terms 6f multiples of the standard deviation Warning = 
Bc+ 3*a^(xVc)"\ Alarm = Bo+ 6*o^ (xSc)"\ It should be noted that a covariance that maybe 
determined as: 



20 



S = s^ (x^)'* where s^ is calculated as noted above. 



As described elsewhere herein, the foregoing techniques are based upon a healthy 
gear characterized as having noise that is stationary and Gaussian in which the noise 
25 approximates a normal distribution. 
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What will now be described are techniques that may be used in detennining an HI 
using the normalized CI values as inputs. In particular, two techniques will be desciibed for 
determining an HI. A first technique may be refeixed to as the non-linear map technique. 
The second technique may be referred to as the hypothesis test method of HI generation. It ' 
should be noted that CI values other than normalized d values may be used in comiection 
v^ith HI determination teclmiques described herein. 

It should be noted that an embodiment may use CI values that are not nomialized in 
connection with the HI determination techniques described hereiiu In this instance, multiple 
torque bands may be used, pne for each CI or group of CIs belonging to different torque 
bands. Additionally, a larger covariance matrix may be used as there may be a larger 
variance causing decrease in separation between classes. 

For any generic type of analysis (gear, bearing, or shaft), a subset of &e 
diagnostics indicators or CIs is selected; The CIs which are best suited to specify the fault 
indication may be developed over time through data analysis. Faults may be calculated at the 
component level and an HI may be calculated for a given component. If there is a 
component fault, then there is a sub-assembly .fault, and therefore a drive train fault 

Following is a description of a non-linear mapping methodology for detennining an 
m. Given a set of component indicators II, 12, 13, „.IN, choose the desired subset of K 
indicators such that K N. For the chosen group of indicators, let WTi define the weight of 
the ith indicator, Wi the warm'ng threshold, and Ai the alarm threshold. Then apply the 
following processing to the set of chosen indicators. 
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Health Indicator Contribution bescrijption 



for XX = 1 JC /* cycle through all K indicators in subset 

If ipoq < Wi /* if less than warning lev^l Wi, assign 0 ♦/ 
Hi contribution = 0 

elseifWi*Ii<Ai 

Hi contribution = l*Wi 

else * 
Hi contribution = 2* Wi 

end 
end 

In the foregoing pseudo-code like description, each indicator or CI is weighted and • 
contributes a portion to the HI determination. Subsequently all the Hi contributions for flie' ' 
selected CIs are summed and may be compared to threshold values for determining one of 
two possible outcomes of "healthy" or "not healthy". 



Consider the following example table of information for a selected subset of 9 CIs 
along with threshold and weight values. It should be noted that in an embodiment, any* one or 
more of the values for weights, warning and alarm values may be modified. 



a 


Value 


Warning 


Alann 


Weight 


HI cbntiibutioii 


No. 




Level 


Level 






n 


3.26 


3.5 


4.0 


1.0 


0.0 


J3 


3.45 


.3.0 


3.5 


1.0 


1.0 


16 


7.5 


6.0 


8.0 


1.4 


1.4 


19 


0.88 


0.5 


0.75 


0.9 


1.8 


114 


4.2 


3-5 


4.5 


1.0 


1.0 


117 


4.7 


3.5 


4.5 


0.9 


1.8 


122 


5.2 


2.0 


4.0 


1.1 


2J2 


123 


4.4 


3.5 


4.5 


1.2 


1.2 


124 


18.9 


10.0 


20.0 


1.0 


1.0 



Using the foregoing example and values, the sum of the HI contributions is 1 1.4. 
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Applying the Health Indicator Contribution technique as set forth in the foregoing pseudo- 
code like description, 12, with a value of 3.26, is below the warning threshold, so the 
contribution to the index is 0. hidicator 13 has a value of 3,45, which contributes a 1 toward 
the index since the weight value is also 1. However, Indicator 16 contributes a 1.4 to the 
indfex because it crosses the warning level (contributing a value of 1 to the index) while being 
weighted by a factor of 1 .4. 

In the foregoing example, if no indicators were in alarm, the smn of HI contributions 
would be zero and if aU indicators were in alarm, the sum would be 19, the worst fault case 
represented by this detector scheme. The HI may be represented as a value of 1 for healthy 
and 0 for not healthy as associated with a component represented by the foregoing C3 values. 

The HI may be determined by dividing 1 1 .4/1 9, the maximum of worst case outcoitae 
to obtain 0.6. This overall health index output ratio can then be compared to another final 
output threshold, where normal components produce His, for example, less than 0.5; values 
between 0.5 and 0.75 represent warning levels, and values over 0.75 r^resent alarm. 

It should be noted that the wei ghts may be determined using a variety of different ' 
techniques. The weights of each CI may be determined using any one or more of a variety of 
techniques. One embodiment may determine weights for the CIs as: 

1 

^ eigen _ values ^of_ihe_ cov ariance _ matrix 

It should be noted that other threshold values may be used in HI determination and 
may vary with each embodiment 
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In one embodiment, using the nonnalized CI described elsewhere herein with the non- 
linear mapping technique, the threshold values may be represented as: Warning «= Bo+ 
3*a^(x*x)"\ Alarm = Bo+ 6*a^(xSc)"*, where Bo may represent a mean or average coe£Eici^t 
as included in the B vector being solved for in the equations described in connection with CI 
normalization. In the foregoing example, the Warning threshold is 3 standard deviations and 
flie Alarm level is 6 standard deviations. It should be noted that other tiu-eshold values may 
be used in and may vary in accordance with each embodiment 

What will now be described is a second technique that may be used in deteramung 
His using CIs, in particular, using nomialized CIs. 

The technique for HI determination may be referred to as Hypothesis testing 
technique for HI determination which minimizes the occurrence of a false alami rate, or 
incoirectly diagnosing the health of a part as being included in the alarm classification when 
in fact the part is not in this particular state. In one embodiment, three classes of health 
indication may be used, for example, normal, warning and alaim classifications with alarm 
being the least "healthy" classification. Other embodiments may use the techniques described 
herein with a different number of classes. As described elsewhere herein, the class of a part 
indicating the health of the part may be determined based on measured vibrations associated 
with the part. Additionally, the technique described herein may use a transforaiation,^ch as 
the whitening transformation to maximize the class distributions or separation of values thus 
decreasing the likelihood or amount of overlap between the classes. In particular, this 
maximization of class separation or distance attempts to minimize the miscla^ification of a 
part. A description of the whitening transformation used in herein in following paragraphs 
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may be found, for example, in "Detection, Estimation and Modulation Theory", Hany h. 
Van Trees, 1968. John Wiley & Sons , New York Library of Congress Catalog Card Number 
67-23331. 



Using the Hypothesis Testing method of HI generation, the HI or classification h(X)" 

of a vector of normalized CI values denoted as X may be detemiihed in which, as discussed 
elsewhere herein in more detail, X may be normalizedUsing the hypothesis testing 
techmque, a determination is made as to which class (nonnal, wanung or alann) X belongs. 
In our instance, there are three classes. However, a first determination using the hypothesis 
testing may be performed using a first class corresponding to normal, and a second class 
corresponding to not noraial. If the determination is normal, then testing may stop. 
Otherwise, if determination is made that the testing results are "not normal", a furthcsror 
second determinatiori using the hypothesis testing may be performed to determine which "not 
normal" class (alarm or warning) X belongs. Thus, the hypothesis testing techmque may be 

performed more than once in accordance with the particular number of classes of an 
embodiment. For three classes, there are two degrees of fi-eedom such that if the sample X is 
not fi-om A or B classes, then it is fi-om Class C. 

X may belong to class tOy or to^, such that: g, (X)^q^ (X) (the notation 5> means 

that if qyPO is greater than q^QC), choose class 2, o>2, or if qy(X) is less than qaCJQ, choose 
class 1, CO,.) In the foregoing, g, is the a posteriori probability of <0/ givenJT, which can be 
computed, using Bayes theorem in which g,= Ppi(X)/pCX), where pCX) is the mixed density 
fiinction. 
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The mixed density ftmction is the probability function for all cases where g i is the 
unconditional probability of "i" given the probability of "i" conditioned on the mixed density 
function. • 

5 Substituting the foregoing representation of eadi ql and q2, since p(X) is common to 

both, now: PxPxiX)oP^p^{X) or as a likelihood fimction as £(^=^^^—^<>-^. The 

likelihood ratio is a quantity in hypothesis test The value P^/P/ is the threshold value, hi 
some instances, it may be easier to calculate the minus log likelihood ratio.* hi this case, the 
decision rule becomes (e.g^ now called the discriminate function): 

10 /j(^) = -hi£(-Sr) = "hij?,(Jr) + hi/7j(A0^hi^ 

jPi 

Assume that the Pi(X)^s are normally distributed with mean or expected values in 
vectors Mr and covariance matrix 2/ . This assumption may be detemiined without loss of 
generality in that, any non-normal distribution can be whitraed, as with the whitening 
transformation described elsewhere herein, with the appropriate power transfonn, or 
1 5 increasing the sample size to the point where the sample size is very large. Given this, the 
decision rule becomes: 
hiX} « - In /(A') 
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Recall tiiat maximization of distance between the two classes is desired to minimize 
the chance of a false alarm or misclassification of a part as broken when it is actually normal. 

A function Z is defined as Z —X-M, (e.g. a shift where Xis the meastired CI data and 
M is the mean CI values for a class), so that: 
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(^) — Z^H'^Z (this distance is the n dimensional distance between two 
distributions). 

Note that S represents the covariance* It may be determined that a particidar Z maximizes the 
distance function, subject to Z^Z = /, the identity matrix. 

Using a standard Lagrange multiplier, y, to find the local extrema (e.g. the maximum) 
a partial derivative is obtained with respect to Z in the following: ' 

' « ' 

where S is the eigenvector of X, 

which may then be set to zero to find the extrema and solving for Z: 
E-'z m pz or zz m xz where X - J/^ ' hi Order that a non-uull Z exits, Xmust be choscn to 

satisfy the determinant: p:^xi\mo - 

Note that X is the eigenvalue of and 2 is the corresponding eigCTvector. S is a 
syinmetric nxn matrix (e.g. a covariance matrix), there are n real eigenvalues (X/.^X,,) and n 
real eigenvectors <pj _ the characteristic equation is: S4E>==4>A, and ^'^<^ = I wher« is an 
« X n matrix consisting of?? eigenvectors and A is a diagonal matrix of eigenvalues (e.g. the 
eigenvector matrix and eigenvalue matrix, respectively). 

y, representing the coordinated shifted value of X, may be represented as: 

having a covariance matrix ofy,2^ = 0^2j^<I> = A where Sx represents the 

covariance of the vector of matrix x . Continuing, the whitening transformation maybe 
defined such that: 

Y = A-^'^<I>^X-(<t>A-^'2)'rX. Sy= A-^^O)^ 2x<l>A-^^= A-^'2^-1/2^^ 
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Thus the transfonnatioii that maximizes that distance between distribution or classes 
is: ' 
A ~Ti^^9^ as shoWn above. 
5 Using this value oif A, define . 

A"^ Si A =1. A"^ S2A = K,and A^ (M2-Mi) = Land ' 
: (sr' 22*')'^ transformed to a diagonal matrix A by A that may be represented as: 

A=A*^[Aa-K-')A'^]-' A=a-K^-)'' ; 

which may be substituted into the discriminate function defined above: 

10 ■ 

Thiis, if the above is less than the threshold, for example. In (P2/P1), then the 
component is a member of the normal or healthy class. Otherwise, the component is 
15 classified as having an HI in the broken class, such as one of alarm or warning. In the latter 
case, another iteration of the hypothesis testing technique described herein may be further 
performed to determine which •'broken" classification, such as alarm or warning in this 
instance, characterizes the health of the component under consideratioiL 

20 In the foregoing technique for h>pothesis testing, values, such as the a posteriori 

probabiUties qi and qa, may be obtained and determined prior to executing the hypo&esis 
testing technique on a particular set of CI normalized values represented as X above. As 
known to those of ordinary skill in the art, Bayes theorem may be used in determining, for 
example, how likely a cause is given that an effect has occurred. In this example, the effect is 
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the particular CI normalized values and it is being determined how likely each particular 
cause, such as a nornial or broken pat, given the particular effects. 

• 

It should be noted that operating characteristics of a system define flje probability of a 
false alarm (PFA) and the probabiUty of detection (PD). The transformation used to 
maximize the distance function optimizes the discrimination between classes. However, the 
threshold value selected given a discriminate function maybe used in detennining the PD and 
PFA. In some embodiments, the cost of a false alarm may be higher than the cost of a missed 
detection. In these instances, the PFA maybe set to define threshold values, and thten accept 
the PD (e.g., a constant false alarm rate (CFAR) type of process). The distance function is a 

normal density function, based on the conditional covariance of the tested vahies under 
consideration. Given that, the PFA may be detennined as: Pp = P(HoHi), which means the 
probability that the sufficient statistic is greater than some threshold is the integral of the 
threshold to infinity of a normal PDF. 

where 

the lower integral limit of 

In this example, the threshold may be the to (P2/P1). This integration is flie 
incomplete gamma function. Conversely, the probability of a detection (PD) is: 



but now 



a = -H^^)/d + and d^ = (M, -M^fLl^M, -Jl/,) 
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Note, the distance function is relative to the condition (e^. Ho or H j) being investigated. 

Referring now to Figure 26, shown is an example of a graphical illustration of the 
5 probability of a false alarm PFA represented by the shaded region A3 which designates fbo 
overlap between the distribution of class HO, denoted by the curve fomied by line Al, and 
class HI, denoted by the curve formed by line A2. 

Referring now to Figure 27, shown is an example of a graphical illustration of the 
probability of an appropriate detection (PD) represented as ai»a A4 as belonging to -cksis 
10 represented by HI as represented by the curve formed by line A2. 

Referring now to Figure 28, shown is a graphical illustration of a relationship in one 
embodiment between the PFA and PD and the threshold value. Note that as the threshold 
increases, the PD increases, but also the PFA increases. If the performaince is not acceptaiUe, 
such as the PFA is too high» an alternative is to increase the dimensionality of the classifier, 
such as by increasing the population sample size, n. Since the variance is related by 
l/sqroot(n), as n increases the variance is decreased and the nomialized distance between the 
distributions will increase. This may characterize the performance of the system. The 
likelihood ratio test used herein is a signal to noise ratio such that the larger the ratio, (e.^., 
the larger the distance between-ttie two distributions), the greater the system pCTformance. 
The process of taking an orthonoimal transformation may be characterized as similar to the of 
a matched filter maximizing the signal to noise ratio. 

Referring now to Figure 29, shown is an example of a graphical illustration of how 
the threshold may vary in accordance with the probability of determining class Ho. 
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It should be noted that false alann rate and detection rate are two factors that may- 
affect selection of particular values, such as thresholds within a particular system. In the 
example embodiment described herein, false alarm rate is a determining factor, for example, 
because of the high cost associated with false alarms and the fact that they may corrode 
confidence when a real fault is detected. It should be noted that other embodim«its and othef- 
applications may have different considerations. Further in this example of the systanof 
Figure 1, certain factors may be considered. An acceptable false alarm rate, for example, 
such as 1 false alarm per 1 00 flight hours, is established. An estimate of the number of 
collection opportunities per flight hours may be determined, such as four data collettions. A 
number of His may be selected for the system, such as approximately 800. A confidence 
level may be selected, such as that there is a 90% probability that a false alarm rate is less 
than 1 per 1 00 flight hours. 

Li this example, it should be noted that each HI is a an independent classification 
event such that the law of total probability may give the system alarm rate using flie 
foregoing: 

System PFA = 1/(100 * 4 * 800) = 3.1250 * 10 

It should also be noted that in the foregoing, when the covariance of two classes is 
approximately the same, or for example, unknown for a class, the logarithm likelihood ratio 
test for classification may be simplified in that the model may be reduced to a linear rather 
than quadratic problem having the following model: 
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If the covaiiance is whitened, the model simplifies further (assuming the appropriate 
transformation is made to the means and measured values). 

What will now be described are techniques that may be used in connection wifli ♦ 
selecting a subset of CIs, such as selection of normalized (Us, for exanaple, under 
consideration for use in determining a particular in. 

If we have a two or more classes (such as alarm, warning and normal classifications), 
feature extraction, or determining which CIs to use in this embodiment, may become a 
problem of picking those CIs or features that maximize class separability. Note tfmt ^ 
separability is not a distance. As described elsewhere herein, an eigenvector matrix 
transformation may be used in maximizing the distance between two fimcfions or distribution 
classes. However, this same technique may not be applicable when some of the infonnation* 
(e,g. dimensionality) is being reduced. For example, in the following test case, three features, 
or CIs, are available, but only two are to be selected and used in determining HI 
classification. The distributions are: 
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When looking at the eigenvalues of the whitening transformation <1.9311,3.0945, 0.4744), 
the maximum distance of the distribution is an axis y (e.g. 2"^ dimension, the distribution was 
whitened and the project dimension (e.g. x, y or z) was plotted), but this axis has the 
minimum separability. Using this as one of the two features will result in higher false alarm 
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rates than another feature. This may identify the importance Of feature 
maximizing the sqjarability. 

The problem of separability may be characterized as a "mixed" problem in that 
differences in means maybe nomialized by diiferent class covariance. If the mean values are 
the same, or the covariance are the same, techmques such as the Bhattachaiyya Distance may 
be used to measure class separability. However, same mean or covariance values may not be 

Ukely and thus such techniques may not be appUcable. Statistical tools developed in 

. ' ■ ■ ■' ' 

discriminant analysis may be used to estimate ciass separabiKty. 

A measure of within class scatter may be represented as the-weighted average of the 
class covariance: = gp,2, . for each class I, where K is the probability of the . 

of the covariance S, for that class. In one embodiment, there may be two classes, such 
healthy or unhealthy. When considering the unhealthy status, for example, when perfonning 
a second round of hypothesis testing described herein, there may be alarm and warning 
classes. 

A measure of between class scatter, Sb, may be represented as the mixture of class 

means: 

Note that Mo represents the mean or expected value of the classes and Mi- Mo is a 
difference or variation from the expected value for the classes under consideration. The 



i occunrence 



as 
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fonnuiation for a criteria for class separability may result in values that are larger when the 
between class scatter is larger^ or when the within class scatter is smaller. A typical critoia 
for this is 7 = diag(S^^S^) , where In general, is not diagonal. One technique takes the 
whitening transformation ofSw where S^A = /, then define the whitening transformation of 
5 Sbas: 

•* - • 

. ' *■ - •; • ■ ■ I " .. 

Now taking the diagonal of the foregoing Sbw gives a better representation of Ac 
class separability of each feature. ' • 

In sumniaiy, CIs may be selected in accordance with the technique described above to 
10 obtain and examine the diagonals of the "whitened" Sb, rq)resented as Sbw. Let X be a 
matrix where rows and colunms represent different ds hiaving a covariance matrix S. An 
embodiment may use normalized CIs and select a portion of ttiese for use. An embodiment 
may also use CIs however, those selected should belong to the same torquie band. 

As described elsewhere herein, let A represent the corresponding eigenvalue matrix 
15 and 4> as the corresponding eigenvector matrix foi* the CI matrix X. Then, A, as described 
elsewhere herein in connection with the whitening transformation, may be repr^esented as: 



where A is the transformation matrix that whitens the covariance E. If Sb is defined as above 
as the between mean covariance of the classes, the whitening matrix A may be used to 
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normalize the differences and give a distance between the mean values of the different 
classes, such that 

Sbw = A'^SbA 

where Sbw represents the *Vhitened» Sb. The diagonals of Sbw may then be sorted m ' 
; descending order in which each diagonal represents an approximation of the size of the 
separation between features or as. Thus, selection ofa subset of*«n** features or CIsfix)m a 
possible **m" maximum Os included in X may be determined by selecting the «n" laigest 
diagonals of the matrix Sbw. In particular, the diagonal entiy 1,1 coiresponds to the first 

column of the covariance matrix and the first a in the vector X, entry 2,2 to the second 
column of the covariance matrix and the second CI in the vector X being considered, and so 
on. 

Once a particular HI is detemiined at a point m time, it may be desired to use 
techniques in connection with trending or predi^g HI values of the component at firtuio 
points in time. Techmques, such as trending, may be used in establishing, for example, when 
maintenance or replacement of a component may be expected. As described elsewhere 
herein, techniques may be used in determining an Hi in accordance with a vector of a values 
having expected CI values included in vector M,- for a given HI classification, i, having a 

covariance matrix S,. One techmque uses a three state Kalman filter for predicting or 
trending fiiture HI values. 

The Kalman filter may be used for various reasons due to the particular factois taken 
into account in the embodiment and uses described herein. It should be noted that other 
systems embodying concepts and techniques described herein may also take into account 
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other noise factors. In one enibodiment, the Kalman filter may be preferr^ in that it provides 
for taking into accoiint the noise of a particular arrangement of components. There may be 
noise corruption, such as indicated, for example, by the covariance matrices described and 
used herein. It may be desirous to filter out such known noise, such as using the Kalman 
filter, providing for smoothing of data values. 

The Kalman filter provides a way to take into account other apriori knowledjge of the 
system described herein. In particular, the health of a component, for qcample, may not 
change quickly with time. The difference between the health of a component at a time t, and 
time t+delta may not be large. This technique may also be used in coimection wifli 
determining future His of a particular part, for example, where the part is old. A part may 
have reached a particular state of relatively bad health, but still a working and functional part 
The techniques described herein may be used with an older part, for exanq>le, as well as a 
newerparL 

In the arrangement Avith the Kalman filter, state reconstruction may be peifoimed 
using the Ricatti equation, as known to those of ordinary skill in the art The technique that 
is described herein uses a three-state Kalman filter of HI, and the first and second derivatives 
thereof with respect to changes in time, denoted, respectively, dt^ and dt^. The Ricatti 
equation in this instance uses a [1x3] vector of time values rather than a single value, for 
example, as may be used in connection with a single state Kalman filter* 

What will now be described are equations and values that may be used in detemiining a 
future value of a particular HI. Let: 
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0 O] 
1 dt 



0 = 



0 0 



F 1 



HI_est 
HI 
HI 



in which: 

o is the power spectral density of the system, 

R is the measurement enxMT, 

P is the covariance, 

^ is the plant noise. 

His the measurement matrix, 

J^is the Kalman gain and 

O is the state transition matrix. 



i 
1 



IS 



H may be characterized as the Jacobian matrix. Since the vahie of a single HI is 
desired, only the first entiy in the H vector is 1 with remaining zeroes. There are n entries m 
the n X 1 vector H for the n state Kabnan filter. Similarly, the X vector above is column 
vector of 3 HI entries in accordance with the three-state Kabnan filter.' The end value being 
determined is the vector X, in this instance which represents a series of HI values, for which 
the first entry, HI_est in the vector X is the one of interest as a projected HI value being 
determined Within the vector X, ii7 represents the first derivative of HI_est and SS 
20 represents the second derivative of HI_est / represents the average amount of time 
between measurements or updating of the HI value. In other words, if dt represents a 
measurement or delta value in units of time between HI deteiminations, and this is perfoimed 
for several instances, ? represents the average of the delta values representing time changes. 
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What will now be presented are equations representing the relationships between the 
above quantities as may be used in determining a value of X(l) for predicting or estimating 
an HI value at a future point in time given a current HI valiie. 

■ 5 ; ■ 

-^lM=<^*^^-I|'-• (Equation Tl) 

"P/k-i =<I>^-.i»-i<I»'"+fi (Equation T2) 

K^P^,H^{HP^,H^-^-R) (Equation 13) 

P^^{J-KH)P^ (Equation T4) 

X^^X^-\-K{HI-HX^) (EquatiOT T5) 

Note that the subscript notation above, for example, such as "tft-l" refeis to 
detennining a value of at a time t conditioned on the measurement at a time of"t-l". 
10 Similarly, "t|t" refers to, for example, detennining an estimate at a time "t" conditioned on a 
measurement of time "t". 

t 

The current HI determined, for example, using other techniques described herein, may 
be input into Equation T5 to obtain a projected value for HI_est, the best estimate of 

15 current HI. To jiroject the expect HI "n" um'ts of time into the future, input the number of 
units of time "dt" into ^ (as described above), and use the state update equation Equation 
Tl) where now Equation Tl becomes: X t+dt|i = * X t|t . This allows the best prediction of 
nicest any numbCT of units of time into the future where HI_est is desired. It should be noted 
that as set forth above, the linear matrix operation such as <E^ X is equivalent to an integration 

20 from t to dt of the state of X, where X represents the vector of HI values set forth above. 

: Different values may be selected for ira'tial conditions in accordance with each 

embodiment. For example, an im'tial value for P representing the covariance may be 
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(1 /mean time value between failures). An embodiment may use any one of a variety of 
different techniques to select an initial value for P. Additionally, since P converges rapidly 
to an appropriate value and the time between data acquisitions is small in comparison to tbe 
mean failure time, selecting a particularly good initial value for P may not be as important as 
5 other imtial conditions. A value for o may be selecting in accordance with qpribn 
information, such as manufacturer's data on the mean time between component parts' 
expected failure time. For example, for a transmission, the mean failure time may be 
approximately 20,000 hours. The spectral deiisity may be set to (1/20,000)*. It should be 
noted that the failure rates may be generaUy characterized as an exponential type of 
10 distribution. The mean time between expected failures is a rate, and the variance is that rate 
to the second power, R may also be deterniined using fl/7rror/ informatioOi such as 
manufacturer's data, for example, an estimated HI variance of manufacturer's data of a 
healthy component. Q may be characterized as the mean time between failures and dt (delta 
change in time between readings). As the value of dt increases, Q increases by the third 
15 power. 

Input data used in the foregoing trending equations may be retrieved from collected 
data, for example, as may be stored in ihc system of Figure 1 . 

JO In determining His, for example, as in connection with the system of Figure 1 for 

particular components. His may be derived using one or more CIs. In calculating Ob, data 
acquisitions may occur by recording observed data values using sensors that monitor different 
components. There may be a need for estimating data used in connection with CI 
calculations, for example, in instances in which there may be too h'tUe or no observed 

5 empirical. For example, in connection with a power train, there may be a need to obtain 
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estimated data, for example, for each bearing, shaft and gear within €ie power train to 
calculate CIs. However, insufScient empirical data may eixist in connection with gear or 
bearing related measurements, such as, for example, those in connection with a^gear or 
bearing related measurements, such as, for example, those in connection wifli a gear or 
S bearing fault due to the rare occurrence of such evmts. In such instances, mean and titireshbld 
values may be deriyed using other techniques. * 

A CI may indicate a level of transmission error, for example, in whidi truismissioii 
eiTor is a meaisure of the change in gear rigidity and spacing. Modeling transmission error 

10 may allow one to gauge C3 seiisitivity and derive tlu-eshold and mean values indict 
gear/bearing failure. This transmission enx>r modeling may be refered to as d>qiamic 
analysis. Whiat will now be described is a technique that may be used to model a geais to 
obtain such estimated values. By modeling eadi gear pair as a damped spring model with ttie 
contact line between the gears, transmission error may be estimated. It should be noted tiiat 

15 this model uses two degrees of jfreedom or movement. Other systems may use other models 
which may be more complex having more degrees of freedom. However, for flie puiposes of 
estimating values, this model has proven accurate in obtaining estimates. Other 
embodiments may use other models in estimating values for use in a system such as that of 
Figure 1. 

20 

Referring now to Figure 30, shown is an example of an illustration of a pair of gears 
for which a model will now be described. A force P at the <:ontact gives linear and torsional 
response to each of the 2 gears for a total of four responses as indicated in Figure 30. TTie 
relative movement dztP is the sum of the 4 responses together with the contact deflection 
25 due to the contact stifiness Sc and the damping coefficient be. This may be represented as: 
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^/sp-^Jcjbp-mpar /fgy-i'jcoqp-Ipa}^ • 



5 



in whidi: 



sp is the linear stiflfiiess of the pinion; , 
j is the square root of -1 ; 

10 CO is the angular rate that may be obtained from the configuration file (e.g., shaft ipm 

* 60 * 2^ to obtain radians per second for the pinion driving the wheel); 

bp is the linear damping coefficient of the pinion; 
mp is the mass of the pinion; 
ip is the radius of fhc pinion; 

15 

kp is the angular efFective stiffiaess of the pinion; 
qp is the angular damping coefficient of the pinion; 
Ip is the angular effective mass of thei pinion; 

20 sw is the linear stiffiiess of the wheel; 

bw is ttie linear damping coefficient of the wheel; 
mw is the mass of the wheel; 
ip is the radius of the pinion; 
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kw is the angular effective stiflBoess of the wheel; 
qw is the angular damping coefBcient of the wheel; 
Iw is the angular effective mass of the wheel; * 
sc is the linear stifQiess of the contact patch where the two gears come into contact; 
be is the linear damping coefficient of the contact patch; 

It should be noted that values for the above-referenced variables on the rights hand 

. I ■ 

I f • ■ ■ . . 

side of EQUATION Gl above, except for P (described below), may be obtained using 

manufacturer's specifications for a particular arrangement used in an embodiment. - . An 

embodiment may include quantities for the above-referenced variables in imits, for exanaple, 

such as stif&ess in units of force/distance (e.g., newtons/meter), mass in kg vaaits, and the 

like. 

The relative movement, d, is the TJE., so from d, the above-referenced equation can 
be solved for P, the tooth force. Deflection is the force (input torque divided by the pinion 
base radius) * the elastic deflection of the shafts, which may be used in ^estimating P 
represented as: 

P - (1/lq) * rp) + (1/sp) + (1/sw) + (l/kwrw)E^^ 
where the variables are as described above in connection with EQUATION Gl. 
Using the above estimate for P with EQUATION Gl, the displacement, such as a vibration 
transmitted through the bearing housing and transmission case (which acts an additional 
transfer function), may be determined. 
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Referring again to Figure 30, shown is an example of an illustration of the gear model 
and the different variables used in connection with EQUATION Gl and G2. Lp may 
represent the longitudinal stiffiiess of the pinion and Lw may represent the longitudinal 
stif&ess of the wheel. It should be noted that these eleinents may not be included in an 
embodiment using the two degrees of freedom modeL 

Bearings may also be modeled to obtain estimates of fault conditions in instances 
where there is little or no empirical data available. With bearings, a periodic impulse is of 
interest The impulse is the result of a bearing rolling over a pit or spall on the inner or outer 
bearing race. The intensity of the impulse on the bearing surface is a iunction of the angle 
relative to the fault, which may be represented as, for example, described in the Stribeck 
equation in a book by TJl. Harris, 1966, Rolling Bearing Analysis. New Yoric: John Wiley p 
148 as: 

q(©)=qot-(j^)[i-«>se)]' equationbi 

where n = 3/2 for ball bearings and 10/9 for rolling elements bearings e < .5, and 6 is less than 
7t/2 in accordance with values specified in this particular text for the different bearings used 
in the above-referenced Stribeck equation represented as in EQUAHON Bl. 

An impulse in a solid surface has an exponential decay constant, which may be taken 
into account, along Avith a periodic system due to rotation of the shaft. The bearing model 
may then be represented as a quantity, "s", which is tiie multiplication of the impulse, "inq>" 
below, the impulse intensity, "q(6)" as may be determined above, the period shaft rotation, , 
which is "cos(e)" below, all convoluted by the exponential decay of the mataial and 
represented as: 
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[imp X g{0)x cos 0]®&cp{T/i) EQUATION B2 

where T is the exponential decay and t is the time. It should be noted that "T" varies with the 
5 material of the solid surface. "exp(T/t)" may be obtained, for example, using a modal 

hammer, to generate the decay response experimentally. An embodiment may also obtain this 
value using other information as may be supplied in accordance with manufacture's 
information. The value of "t" may be a vector of times starting with the first time sample md 
extending to thfe end of the simulation. T is generally small, so the expression "exp(T/t)" 
10 approaches zero rapidly even using a high sampling rate, \ . ^ ' ■ 

"imp** is the impulse train that may be represented as the shaft rate ^ bearing frequncy 
ratio * sampling rate for the simulation period, 

15 "s" is the simulated signal that may be used in deteraiining a spectrum, "S", where "S 

= fft(s)", the Fourier transform of s into the frequency domain from the time domain. As 
described in more detail in following paragraphs, in determining a CI in connection with the 
bearing model signal "s" having spectrum "S", for example, the Power Spectral Density of S 
at a bearing passing frequency may be used as a CI. Additionally, for example, other CI 

20 values may be obtained, such as in connection with the CI algorithm comparing the spectrum 
"S" to those associated with transmission error in connection with a nomial distribution using 
the PDF/ CDF CI algorithms that may be generally described as hypothesis testing techniques 
providing a measure of difference with regard whether the spectrum is normally distributed. 
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It should be noted that, as described elsewhere herein in connection with gear models 
values may be used in the foregoing equations in connection with simulating various fauK 
conditions and severity levels. The particular values maybe detennined in accordance wilh 
what small amount of observed data or manufacturer's data may be available. For example, 
in accordance with observed values, an impulse value of 0.02 for the impulse, "imp", may • 
correspond to a fairly severe fault condition. Values ranging fiom'O.OOl to 0.03, for exanq>le, 
may be used to delimit the range of "imp" values used in simulations. 

Following is an example of estimated data using the foregoing equations for a bearing 
having the followmg configurations: 

Rpm = 287.1 
Roller diameter B.2S 
Pitch diameter «= 1.4171 
Contact ai^e = 0 
Number of elements = 10 
Inner race fault 

Refraxing now to Figure 31, shown is an example of a graphical representation of ^ 
signal for the foregoing configuration when there is some type of bearing fault as estimated 
using the foregoing equations EQUATION B 1 and B2. Figure 32 represents the estimated 
spectrum "S" as may be determined using EQUATION B2 above. 
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It should be noted that for bearings, there may be three typeis of faults/ for example^ 
estimated using the foregoing equations. There may be an inner race fault, an outer race fault 
or a roller element fault Localized bearing faults induce an excitation which can be modeled 
as an impulse train, expressed as imp . in the above equation. This impulse "imp" corresponds 
to the passing of the rolling elements of the faiilt Assuming b. constant, inner ring rotatioii 
speed, the impulse train is periodic and the periodicity depends oii the fault locadaiL 

For outer race faults, the bearing frequency ratio, f a, or may ^ represented as: 

/.^=Y[l-|^<x>s(a)j(/^--/^) EQUATIONB3 



where: 

"db" represents the roller diameter. 

It 

"dm" represents the pitch diameter, 
15 "a" * frequency,>i. 

"fir " is the rotation frequency of the inner race (e.g. shaft rate), and 
'%r " the rotation frequency of the outer race (if fixed == 0), 



20 For inner race faults, the bearing frequency ratio, f a. ir may be represented as: 



= Y^l+|*-cos(a)j(A -far) EQUA'nONB4 



Replacing a with 2nf^^ the time response is fft). This substitution may be performed 
as the initial value of a is based on an angle and not a function of tiaie. In a simulation, there 
25 is a time dependent response as expressed using f(t). 
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The radial load applied to the bearing is not constant and results in a load distribution, 
which is a function of angular position. If the defect is on the outer race, the amplitude of the 
impulse is constant because the fault location is not time varying. For an inner race fault, the 
amph'tude with respect to angular position. The function is: 




EQUATION B5 



9(0 = q(27f(0)). This quantity q(t) is amphtude at a. particular time, or q(theta) representing 
the amplitude at a particular angle. Amplltuide modulation takes into account the distance 
from the fault to the sensor. For outer race fault, the quantity cos (6) is constant (1), for inner 
race fault, it is the cosine function, noted as "cos (6)" in the above equation. 

For a linear system, the vibrations at a given frrvjuency naay be specified by the 
amplitude and phase of the response and the time constant of the exponential decay. As the 
angle, 6 above, changes, the impulse response, h(t), and the transfer function H(f) also 
change due to the changing transmission path and angle of the applied impulse. It is assumed 
that the exponential decays is independent of the angle 0, so that the response measured at a 
transducer due to an impluse applied to the bearing at the location 9 is characterized by an 
amph'tude which is a function of G. 

The impulse response function h(t) and the transfer function H(f) may be replaced by 
a function a(0) giving the amplitude and sign of the transfer function H(f) 
at each angle theta and by the exponential decay of a unit impulse, (e(t)). For an inner race 
defect, rotating at the shaft frequency fs, the instantanous amplitude of the transfer function 
between the defect and the transducer as a function of time, a(t) may be obtained by 
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substituting 27i*fs*t for theta. Note that a(t) is periodic. At 6 = 0 relative to the defect and 
transducer, a(t) has its maximum value. At 6 = it, a(t) should be a minimum because ihc 
distance fonn the defect to the transducer is a minimum. Additiondly, the sign is negative 
because the impulse is in the opposite direction. Because of these properties, the cos(t) may 
5 be used for the function a(t). 

The impulse train is exponentially decaying. The decay of a unit impulse can be 
defined bjr: , , ' 

10 e(0 = exp(-^) EQUATION B6 

for f> 0, where Te is the time constant of decay* 

The bearing fault model is then: 

v(/) = [//?/p(0^(/M/)]*<0 EQUATION B7 

where: 

imp(t), which is the impulse over a time t, = 2ic * shaft rate * tune * bearing 
frequency ratio, as may be determined using EQUATIONS B3 and B4 above; 

a(t) is the cos(e) for an inner race, which is 1 for an outer race, where cos(6) — 0, 
where 0 is time varying; 

and q(t) and e(t) are as described above. 

25 An embodiment may include a signal associated at the sensor for gear and bearing noise 

combined from the bearing and the gear model may be represented as: 
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s(t) = [d(t)f(t)git)aiO]*e(t)*h(t) EQUATION B8 

where; 

hCO is the frequency response of the gear case, as may be detenmned, for sample, 
using an estimate produced with linear predictive coding (LPC) techniques or wifli a modal 
hammer analysis; 

<f^y is the signal associated with gear/shaft TJB. as may bfe determined using the gear 
model EQUATION Gl; 

and other variables are as described elsewhere herein. 

The frequency spectrum of signals representing a combined bearing and gear model 
fipm EQUATION B8 may be represented as: 

S(f)^\p{f)*Fif)*Qif)*A(J)]E(J)Hif) EQUATIONB9 

As described elsewhere herein, healthy data, such as may be obtained nfiwi£ 
manufacturer's information, may be used in determining different values, such as fliose in 
connection with stiffiiesses for gear simulation, amplitude and exponential decay forbearing 
faults. In terms of generating fault data, since these systems are linear, the following may be 
defined: 

• For gear faults indicative of a cracl^ a reduction in the stiffiiess for a toofli (e.g. 50 and 20 
percent of normal) may be used in estimating median and high fault values. Additionally, 
these values may be varied, for example, using the Monte Carlo simulation to quantify 
variance. 

• Fbr shaft misalignment, shaft alignment within the model may be varied to estimate 
mean fault values 
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• For gear spalling faults, the "size" of an impulse may be detennined through trial and 
error, and by comparing simulation values with any limited observed fault data previously 
collected^ t 

• For bearing fault models, which are spalling faults, fhe size of an impube, indicative of a 
fault, with known bearing faults, may be detennined similarly as v^tb ^ear spalling faults 

Sensitivity analysis may be performed, for example using range of different inpvt 

values for the different parameters, to provide for increasing the effectiyen^ of fault 

.**-■'». 

detection teclmiques, for example, as described and used herein. For example, an * 
embodiment may be better able to simulate a family of bearing faults to iaQor a particular d 
algorithm to be sensitive to that particular fault 

Using the foregoing, the modulated transmission error of a gear mesh, for sample, 
which is a signal may be simulated or estimated. This signal may ^bsequently be processed 
IS . usinganyoneormoreofavariety of CI algorithms such .&at estimates for the mean and 

threshold values can then be derived for fault conditions. (It is assumed that the stiffiiess and 
torque are known apriori). Parameter values used in the above equations corresponding to a 
healthy gear, for example, as may be specified using manufacturer's data, may be modified to 
estimate parametier values in connection with different types of faiiksbeiiigsm By 
20 modifying these parameter values, different output values may be determined corresponding 
to different fault conditions. 

For example, known values for stifBiess, masses, and fhe like used in EQUATION Gl 
may be varied. A cracked gear tooth may be simulated by making the stijGfiiess time varying. 
25 The contact pitch may be varied with time in simulating a shaft alignment fault A modulated 
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input pulse on d may be used in simulating a spall on a gear tooth. Different parameter 
values may be used in connection with specifying different degrees of fault severity, such as 
alarm levels and warning levels. A particular parameter value, such as a tooth stiflBiess of 
70% of flie nomial manufacturer's specified stififiiess, may be used in simulating waraing 
levels. A value of20%ofthenoimal manufacturer's specified stiffiiess may be used in 
simulating alarm levels. The particular values may be determined* in accordance wifli 
comparing calculated values wifli the characteristics of real CI data on any few real faults 
collected. , 

While the invention has been disclosed in connection with the preferred embodiments 
shown and described in detail, various modifications and improvements thereon will become 
readily ^parent to those skilled in the art Accordingly, the spirit and scope of the present 
invention is to be liniited only by the following claims. 
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What is Claimed isr 

1 . A method executed in a computer system for detemuning a health indicator 
associated wifli a component comprising: 
5 determining a plitrality of health classifications; 

determining at least one condition indicator quantifying a 6hiaracteristic of said 
conq>onen^ 

determining a probabihty associated with each of said health classifications^ said 
probabflity being an estimation that said component is of a particular health* classification 
10 given said at least one condition indicatOEr; and 

determining, for a given set of observed values, which of said plurality of health 
classifications is associated with said component usinjg said probabilities associated with said 
health classifications, 

15 2. The method of Qaim 1, wherein said detemiining at least one condition indicator, 

said determining a probability, and said determining which of said health classifications is 
associated with said component, are perfomied upon two of smd.pliiraUty of he^ 
classifications, axid 

wherein, if a determination is made that a first of said health classifications is not 
20 associatedLW^^^ determining at least one condition indicator, 

said determining a probability, and said determining which of said health classifications is 
associated with said component using two other health classifications, 

25 
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3. The method of Claim 2, further comprising: 

storing previously observed and calculated data associated with said component; and 
using a portion of said previously observed and calculated data to determine said 
probability associated with said each heallh classificatioiL 

4. The method of Qaim 2, further comprising: 

obtaining data estimates using models of said component; and 
using said data estimates to determine said probability associated with said each 
healtti classification. 

5. The method of Claim 1, further comprising: 

determining a threshold using a ratio of a portion of said probabilities; 
calculating a current quantity using said observed values; and 

comparing said current quantity to said ratio to determine which of said health 
classifications is associated with said componenl 
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. 6. The method of Claim 5, wherein said calculating said current quantity further 
comprised: 

detenniiiing a transformation matrix maximizing a distance be^ 
plurality of health classifications; 

determining a covaiiance matrix for each of said health classifications using said at 
least one condition indicator and said observed data; and V 

determining said current quantity in accordance with said covariance matrices^ said 
transformation matrix and said observed data* 

7. The method of Claim 6, wherein, for two of said plurality of health classifications, 
said current quantity^ hpO, may be represented as: 

MJ0=|y'"A-'3^-[(-K-'L)']r+[~rK-'i-|ln|K|-^^^ 
wherein* 

X represents said at least one conditional indicator forming a vector of individual 
conditional indicator values; 

Y represents A"^'^ a)^X ; 

A rq)resents (I - K'*)'\ diagonal matrix of eigenvalues of X using a characteristic 
equation of S4> = «>A,^ ^ = identity matrix; 

„ represents an nxn nciatrix of eigenvectors^ 4>1 «^^»iiH>f ^-^^^ 

equation; 

I represents A EiA; 

K represents A^SaA; 

L represents A'^^CMZ - Ml); 
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P2 represents an aposteriori probability of a first of said two healtb classifications 
given : - ■ 

PI represents an aposteriori probability of a second of said two health classifications 
given Si 

5 A rq^resents said transformation matrix maxi^^ 

of said two health classijQcations represented as A"''^***' ; 

51 represents a covariance matrix of said first health classificati^^ 
expected values Ml; and 

'■■ ' ■ ■ ' 

52 represents a covariance matrix ofsaid second health classification having a ve^ 
10 of expected values M2. 

8. The method of Claim 1, wherein said at least one condition indicator indicates a 
physical state of a portion of said component 



82 



wo 02/095633 PCTAIS02A16380 
9. A method for determining a health status of a component compiistng: 
selecting a plurality of condition indicators having a value and each having a 

corresponding weighting factor, and at least one threshold value definiiig at least two 

classifications; 

5 deterimning a contribution to a health indicator for eadi of said condition indi<^^ 

wherein said determining further comprises, for each of said pi 

determining which of said at least two classifications said value of said eadb. 
indicator belongs; and 

determining said contribution to said health indicator by said' eadi cdnditioai 
10 * indicator in accordance with a selected one of said at least two classifications and said 
weighted value; and ^ 

determining said health indicator in accordance with all contributions by each of said 

tt ■ . ■ 

condition indicator values. 

15 10- The method of Claim 9, wherein three classifications are associated with each of, 

said condition indicators formed by two threshold values, and the method fiirther comprises: 
determining that said contribution to said health indicator for a first of said condition 
indicators is zero if said value of said condition indicator is in said fij^ classification; 

determining that said contribution to said health indicator is a first multiple of said 

20 weighting factor if saidjvaluejof-said condition indi 

determining that said contribution to said heal& indicator is a second multiple of said 
weigjiting factor if said value of said condition indicator is in said third classification. 

11. The method of Claim 1 0, where the two threshold values are alami lev^ aiul 
25 warning level. 
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12. A computer program product for determining a health indicator associated with a 
component comprising machine executable code for 

determining a plurality of health classifications; 

5 ^determining at least one condition indicator quanting a characteristic of said 

component 

determining a probabihty associated with each of said health classifications, said 
probabiUty being an estimation that said component is of a particular health classification 
given said at least one condition indicator; and 
J deteimining, for a given set of observed values, which of said plurality of healtfa 

classifications is assodated with said component using said probabilities associated with said 
health classifications. 

13. The computer program product of Qaim 12, wherein machine executable code 
for detennining at least one condition indicator, determining a probabiKty, and determining 
which of said health classifications is associated with said component, arc performed upon 
two of said plurality of health classifications, and the computer program product further 

includes machine executable code, wherein, if a determination is made that a first of said 
health classifications is not associated with said component, for repeatedly deteimining at 
least one condition indicator, repeatedly determining a probability, and repeatedly 
detennining which of said health classifications is associated with said component using two 
other health classifications. 

14. The computer program product of Claim 13, fiirther comprising machine 
executable code for: 

storing previously observed and calculated data associated with said component; and 
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using a portion of said previously observed and calculated data to determine said 
probability associated with said each health classification. 



15. The computer program product of Claim 1 3, further comprising machine 
S executable code for 

obtaining data estimates using models of said component; and 
using said data estimates to determine said probability associ^ed with said each ; 
health classification. 

10 16. The computer program product of Claim 12, further comprising machine 

executable code for ; 
determining a threshold using aratio of aportion of said probabilities; 
calculating a current quantity usrag said observed values; and 
comparing said current quantity to said ratio to determine which of said health 

15 classifications is associated with said component 

17. The computer program product of Claim 16, wherein said machine executable 
code for calculating said cunrent quantity further comprises machi^ 

determining a transformation matrix maximizing a distance between two of said 
20 plurality-of4aealfli classifications; - - — ^ 

determininig a covariance matrix for each of said health classifications using said at 
least one condition indicator and said observed data; and 

determining said current quantity in accordance with said covariance matrices, said 
transfoimation matrix and said observed data. 

25 
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1 8. The computCT program product of Claim 17, wherein, for two of said plurality of 
health classifications, said current quantity, hpO, may be represented as: 




wherdn, . 



X rqjresents said at least one conditional indicator foiming a vector of individual 
conditional indicator values; 

Y represents A'^ *TX ; 

A represents (I - K"*)"', diagonal matrix of eigenvahies of X using a characteristic 
equation of E4> = <i>A, = identity matrix; 

4> represents an nxn matrix of eigenvectors, fl>l.. <Pn, of X using said characteristic 
equation; 

I represent^' A^SiA; 

K represents A^SzA; 

L represents A^(M2 - Ml); 

P2 represents an aposteriori probability of a first of said two health classifications 
given X; 

PI represents an aposteriori probabiKty of a second of said two heallfa classificatiOTis 
given X; 

A represents said transformation matrix maximizing a distance between distributions 
of said two healtb classifications represented as 

Si represents a covariance matrix of said first health classification having a vector of 
expected values Ml ; and 
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22 represents a covariance matrix of said second health, classification having a vector 
of expected values M2. 

19. The computer program product of Claim 12, wherein said at least one condition 
indicator indicates a physical state of a portion of said component 
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20. A computer program product for d etermining a health status of a component 
comprising machine executable code fon 

selecting a plurality of condition indicators having a value and each having a 
corresponding weighting factor, and at least one threshold value defining at least two 
classifications; * 
determining a contribution to a health indicator for each of said condition indicators^ 
v/herein said machine executable code for determining further comprises, for each of said 
plurality of indicators, machine executable code for 

detemiining which of said at least two classifications said value of said each 
indicator belongs; and 

determining said contribution to said health indicator by said each condition 
indicator in accordance with a selected one of said at least two classifications and said 
weighted value; and 

determining said health indicator in accordance with all contributions by each of said 
condition indicator values. 

21. The computer program product of Claim 20, wherein three classifications are 
associated with each of said condition indicators formed by two threshold values, and the 
computer program product fiirther concq^rises machine executable code fon 

dete rminin g that said contribution to said health indicator for a first of said condition 
indicators is zero if said value of said condition indicator is in said first classification; 

determining that said contribution to said health indicator is a first multiple of said 
weighting factor if said value of said condition indicator is in said second classification; and 

determining that said contribution to said health indicator is a second multiple of said 
weighting factor if said value of said condition indicator is in said third classification. 
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22. The compute program product ofQaim 21, where the two threshold valu« are 
alann level and wanung leveL ' 
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23. A method executed in a computer system for determiniiig an health mdicator of a 
component at a subsequent time comprising: 

detemiining a first health indicator of said component at a time^ n, in accordance witli 
at least one corresponding condition indicator and * 

using a three state Kalman filter to determine a second health indicator of said 
component at a time subsequent to time n. 

24. The method of Qaim 23, fiulher comprising: 

estimating said first health indicator using a hypothesis determination technique. 

25. Themethodof Claim 24, \vherein the health of a component does not vaxy by a 
large amotmt with respect to a change in time. 

26. Hie method of Claim 24, wherein said at least one condition indicator is a 
nomialized condition indicator. 
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27. The me&odofClaim 23, wherein 
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in which: - 
a is a power spectral denaly, 
R is a measumnent enxtr, 
P is a covariance* 
j2isap]aiitnoi8p, 
^ is a measurement matrix, 
J^isaKahnan gfiin and 
(I> is state transition matrix, and 

K = P^H^iHP^H^ +/?) 

P^=iI-KH)P^^ 

^^^X^+K{HI-HX^) 
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(Equation Tl) 
(Equation T2) 
(Equation . T3) 
(Equatitm T4) 
(Equation T5) 



and the mediod comprising: . 

using flie above Equations; Tl ftrbugh TS to detennine a value Hi__est whidi is a 
estimate of a value HI at a time in the future. 
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28. A computer program product for determining an health indicator of a component 
at a subsequent time comprising machine executable code fixn 

detemiining a first health indicator of said componrat at a time» n, in accordance with 
at least one corresponding condition indicator; and 

using a three state Kalman filter to determine a second health indicator of said 
conotponent at a time subsequent to time n. , 

29. The computer program product of Claim 28, further comprigfng machine 
executable code for 

estimating said first health indicator using a h>pothesis detemiinatian technique. 

30. Thecomputerprogramproduct of Claim 29, wherein the healtib of a conq>one^ 
does not vary by a large amount with respect to a dbange IB tinie. 

31. Thecomputerprogram productofClaim 29, wherein said at least one condition 
indicator is a normalized condition indicator. 
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32. The computer program product of Oaim 28, wherein 



H = [l 0 O] 

1 dt 



0 0 



fi = 2CT»r 



? 1 



t 
1 



m 
in 



in which: 

<T is a power spectral density, 
R is a measurement eoor, 
P is a covariancep 

JET is a measurement matrix, 

J^is a Kahnan gain and ■ 

<I> is state transition matrix, and 



(Equadon Tl). 

<Equatioa T2) 
(EquatidQ T3) 
{Equation T4) 

^Bquation TS) 



and the computer program product coniprises machine exiscutable code for: 
detominijig, using above Equations Tl through T5, to determine a value HSjest 
which is a best estimate of a value HI at a time in the future. 
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33. A method executed in a computer system for ranking condition indicators used in 
determining a health indicator for a component comprising: 

determining a first set of a plurality of said condition indicators; 

determining a covariance matrix corresponding to said plurality of condition 
indicators; 

detemiining a transfonnation matrix that whitens the covariance matm 

using said whitening matrix to detennine differences between said first plurality of 

condition indicators and expected values for said condition indicators belonging to a. health 

class» each health class having a corresponding health indicator; and 

selecting a portion of said plurality of condition indicators in ac(x>rdance 

condition indicators have the smallest of said differences. 

34. Themethodof Claim 33, further comprising: 

sorting said differences in descending order, each of said differences having a 
corresponding condition indicator. 

35. The method of Qaim 33, wherein said first plurality of condition indicators 
coirespond to an observed data acquisition. 
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36. The method of Claim 33, further comprising: 
determining a measure of between class scatter, Sb, represented as: 

where MO is an expected value of all L classes. Mi represents an expected value of a 
5 particular class, and Pi is the probability of a class L 

37. The method of Claim 36, further comprising: 

detennining a whitening transformation of Sb as Sbw represented as: 
Stw^A^ SiAt for the whitening transformation matrix A. 

>. - ■ ' \ ■ 

38. The method of Claim 37, wherein said whitening transfbmiation matrix A is: 
^-I'^^T ^ corresponding eigenvalue matrix A and a corresponding eigenvector 

matrix 4>. 
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39* A computer program product for ranking conditiOD indicators used in determixmig 

ahealth indicator for. a component comprising machine executable code foR ^ . ^ 
determining a first set of a plurality of said condition indicaton; ^ 
determining a covariaiice matrix correqp|onding to said pli^^ 

indicatois; 

determining a transformation matrix that whitens the coyanance matrix; 

using; said whitening matrix to determine differences between said first plur^ty of 
condition indicators and expected values for said condition indicators belonging to a beatlli 
class, each health class having a corresponding health indicator; and i / , 

selecting a portion of said plurality of condition indicators in accpFdance with fbMC 
condition indicators have the smallest of said difTerences. 

40. The coniputer program product ofClaim 39, finihcr comprising machiM 
executable code fijn 

sorting said differences in descending order, each of said differences ha^^ 
coiresponding condition indicator. 

41. The computer program product of Claim 39, wherein said first plurality of . 
condition indicators correspond to an observed data acquisition. 

42. The computer program product of Claim 39, fijrfber comprisu^ machine 
executable code for 

determining a measure of betweejn class scatter, Sb, represented as: 
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where MO is an expected value of all L classes. Mi represents an expected value of a 
particular class, and Pi is the probability of a class i. 

43- The computer program product of Claim 42, further comprising machine 
executable code fim 

determining a whitening transformation of Sb a$ Sbw represented as: 
-^^iSJfr^, for the whitening transfomiation 

* 'I 
44. The computer program product of Claim 43, wherein said whitening 

transformation matrix A is: 

A"'^fl>^ with a conre^onding eigenvalue matrix A and a corresponding eigenvector 

matrix 
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45. A method for estimating a conditional indicator value associated widi a gear pair 
comprising: 

5 modeling said gear pair as a damped spring model having a contact line between said 



deterniining a force, P» at a point of contact along said contact liM causii^ 
torsional response to each of said two gears in said gear pair; 

detemiining a relative movement, d, of said gear pair, in accordance wifli said fi»ce, 
10 P, as a smn of four responses and a contact deflection, said relative movement d representing 
a gear model hiaving two degrees of freedom; and 

using said relative movement, d, in determining said conditional indicator value for 
transmission error astociated with said gear pair. 

15 46. The method of Oaim 45, further comprising: 

varying parameter values in accordance with simulating different &ult conditions; and 
determirung mean and tlu-eshold values for said difiFer«nt fault conditions. 
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47- A method executed in a computer system for estimating a conditional indicator 
value associated with a bearing comprising: 

deteraiining a bearing frequency ratio for said bearing; , 
determining a periodic impulse in accordance with said bearing frequency ratio; 
detennining an intensity of an impulse on a bearing surface as a iunctipn of an an^e * 
relative to a bearing fault ^ 
determining a decay of a imit impulse; and 
determining a movement of said bearing^ 

determining a conditional indicator value associated with said bearing in acck>ixlance 
10 with said movement 

48. The method of Claim 47, wherein said bearing frequency ratio corresponds to one 
of: an mner race frequency and an outer bearing race frequency. 

49. The method of Claim 47, finfher coiAprising: 

varying parameter values in detennining conditional indicator values associated with 
different health classifications of said bearing. 

50. The method of Claim 49, fiirther comprising: 
estimating values associated with a gear noise; 
estimating values associated with a bearing noise; and 

estimating a condition indicator of a gear and beating noise as a combined signal • 
using said values associated with said gear noise and said bearing noise. 

51. The method of Claim 50, further comprising: 
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detemiining said combined signal s(t) as: 

^(0 = [rf(0/(Og(0fl(O]* «(0 ♦ A(0 

■ • 

. wbere: ■ " . • * ' 

is a frequency response of a gear case; 
d(t) is a signal associated with gear and shaft transmission enoi; 
f(t) is a bearing frequency ratio; 
q(t) is an amplitude at a particular time t, 

• • • 

a(t) is a cosine for one of an inner and outer race condition at a locatibn theta at whicii 
an impulse is applied to said bearing azid 
e(t) is a decay rate of unit impulse; 

S2« The method of Claim 51, wherein h(t) is determined using one of: a linear 
predictive coding (LPC) technique and modal hanimer anal 

53. The method of Claim 51, iurtfaer comprising: 

detennining a frequency spectrum of signals of s(t) as S(f) in the frequency domain 
by fomung: 

in the frequency domain. 

54. A computer progran product for estimating a conditional indicator value 
associated Avith a gear pair comprising machine executable code fixn 



101 



wo 02/095633 PCTAJS02/16380 

modeling said gear pair as a damped spring model having a contact line between said 

gears; 

determining a force, P, at a point of contact along said contact line causing linear and 
torsional response to each of said two gears in said gear pair; 
5 determining a relative movement, d, of said gear pair, in accordance with said force, ' 

P, as a sum of four responses and a contact deflection, said relative movement d representing 
a gear model having two degrees of fieedom; and 

using said relative movement, d, in determining said conditional indicator value for 
transmission enx>r associated with said gear pair. 

10 

55. The computer program product ofClaim 54, further comprising mach^ 
executable code fan 

varying parailmeter values in accordance with simulating difTerent fault conditions; and 
determining mean and threshold values for said different fault conditions. 

15 

. 56. A computer program product for estimating a conditional indicator value 
associated with a bearing comprising machine executable code for: 
determining a bearing frequency ratio for said bearmg; 

detennining a periodic impulse in accordance with said bearing frequency ratio; 
20 determining an intensity of an impulse on a bearing surface as a function of an angle 

relative to a bearing faul^ 

determining a decay of a imit impulse; and 
detennining a movement of said bearing; 

determining a conditional indicator value associated with said bearing in accordance 
25 with said movement 
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57. The computerprogram product ofQaim 56, wherein said beari^ 
corresponds to one 6£ an inner race fi'equency and an outer hearing race firequency. 

5 58. The computer program product of Claim 56, further comprising machine 

executable code fan * 

varying parameter values in determining conditional indicator values associated witii 
different health classifications of said bearing. 

10 59. Thecomputerprogramproduct of Qaim 58, further compnsingmacl^^ 

. executable code for. ■ . ^ . 

estimating values associated with a gear nois^ 
estimating values associated with a bearing nois^ and 

estimating a condition, indicator of a gear and bearing noise as a combined signal 
15 using said values associated with said gear noise and said bearing noise. 

60. The computer program product of Qaim 59, fuiiher comprising machine 
executable code fim 

determining said combined signal s(t) as: 

20 

where: 

h(i) is a fi-equency response of a gear case; 
d(i) is a signal associated with gear and shaft transmission enor; 
25 f(t) is a bearing fi^uency ratio; 
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q(t) is an amplitude at a particular time t, 
a(t) is a cosine for one of an inner and outer race condition at a location theta at which 
an impulse is applied to said bearing; and , 
e(t) is a decay rate of unit impulse; 

61 . The computer program product of Qaim 60, wherein li(t) is detemiined using one 
of: a linear predictive coding technique and modal hammer analysis. 

' I * ■ 

62. The computer program of Claim 60, further comprising madiilie executable code 

determining a frequency spectrum of signals of s(t) as S(f) in the frequency domain 
represented as 

in the fi^uency domam. 
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63. A method for nonsalizihg a set of at least one observed condition indicator 
comprising: 

detennining a plurdity of conditional indicators and at least one associated factor in 
accordance with previous data acquisitions; 

determining a mean and at least one model coeiffident coxiesponding to said at least 
one associated factox; and 

adjusting said set of observed condition indicators in accordance with model 
coefficients aiid said at least one associated factor producing a noiinalized set. of condition 
indicators. 

64. The method of Claim 63, wherein said relevant £ulor5 in^ 
torque^ airspeed, and rotor speed. 

65. The method of Claim 64, wherein said relevant factors include torque, airspeed 
and rotor speed (Nr), and the method further con^rising: 

solving for B using the equation fomi y = B x for the following: 



B = [5o Bi..Jm ] and x^ 



1 ti Nk^ Air^>eed\ 
1 U.. iVx... Airspeed... 
1 Nm Airspeedri 



wherein the vector y includes condition indicator values and the matrix x includes relevant 
factors corresponding to one of said condition indicators included in y, and B includes 
coefficients. 
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66. The method of Claim 65, wherein a normalized CI is determined at a later point in 
titne by determining: 

CInonnalized = Clobs - (B ♦ x) 

r 

where Clobs is an actual observed value corresponding to a nomlalized CI vahie 
Clhormalized* . 

* ■ * ■ 

67. The method ofClaim 66, furflier comprising: 
storing said values for B and X in a database: and 

retrieving said values for B and x in connection with determining Claormalized. 

68. A computer program product fox normalizing a set of at least one observed 
condition indicator comprising machine executable code for: 

determining a plurality of conditional indicators and at least one associated factor in 
accordance wiHi previous data acqidsitions; 

determining a mean and at least one model coefficient corresponding to said at least 
one associated factor; and 

adjusting said set of observed condition indicators in accordance with model 
coefficients and said at least one associated factor producing a normalized set of condition 
indicators, _ 

69. The computer program product of Claim 68, wherein said relevant factors include 
at least one of: torque, airspeed, and rotor speed. 
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70. The computa^ program product of Claim 69, wherein said relevant factors include 

torque, airspeed and rotor speed (Nr), and the computer program product further comprising 

machine executable code fon 

solving for B using Ihe equation fonn y<= B x for the following: 



y = 



B = [Bo Bu.JBm ] and x = 



1 t% N&i Airspisedl 
1 Nm... Airspeed. 

1 Nm Ap'speedn 



wherein the vector y includes condition indicator values and the matrix x includes rekvaot 
factors corresponding to one of said condition uCidicators included in y, and B includes 
coefficients. 

10 71. The computer program product of Qaim 70, wherein a normal^^ 

determined at a later point in time, and the computer program product furthtf compriang 
machine executable code for determining: 

Onormalized « Clobs - (B * x) 
where Clobs is an actual observed value corresponding to a normalised CI value 

15 CInonnalized. 



20 



72, The computer program product of Claim 71 , further comprising machine 
executable code fer: 

storing said values for B and x in a database; and 

retrieving said values for B and x in connection with determining CInormalized. 
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73, A method executed in a computer system for detemnning a con^ 
characteristic of a component comprising: 

detennining a distribution of observed data associated with said component ' 

measuring a difference between said distribution and a norma] distributioii; and 

detennining said condition indicator using said difference^ 

74. The method of Claim 73, iurther comprising: 

determining whether said distribution of observed data is normally disbibut^ using 
said difference using at least normahty test that is one of: chi-square goodness of fit test, 
Kohnogoiov-Smimof goodness of fit test, Lilliefors test of normality and Jaique-Bera test of 
noimaliiy. 

75. The method of Claim 74, wherein said normal distribution is one nf « n nrtn^l 
cumulative distribution fimction and a normal probability distribution fimction in accoidance 
with said at least one normality test 

76. The method of Claim 75, wherein said distribution of observed data associated 
with a component approximates one of: a Gaussian distribution if said conq)onent is healthy 
and a non-Gaussian distribution otherwise. 
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77. The method of Claim 74, further compriang: 

determining a number of differences between said observed data and expected data, 
said expected data being represented by said normal distributioi^ and 
determining a sum using said difTerences; and 

if said number of differences is greater than a critical value, determining that said 
observed data is not nomially distributed, said critical value being a ^ireshold. 

78. The method of Qaim 77, further compzising: 
deteimining a score bdng a maximum deviation fix>m said oriti 

condition indicator being said score. 

79. The me&od of Claim 78, wherein sexisitivity of said condition indicator mcreases 
as a number of observed data values increases. * 

80. The method of Claim 73, wherein said normal distribution qyproximates a 
distribution of expected values. 
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8 1 . A computer program product for determining a condition indicator about a 
characteristic of a component comprising machine executable code fim 

determining a distribution of observed data associated wifli said component; 
measuring a difference between said distribution and a normal distribiiticxn; aid *. 
detemiining said condition indicator using said difTerencGi 

82. The computer program product of Claim 81, further compxisiog: 

machine executable code for determining whether said distribution of observed data is 
normally distributed using said difference using at least normality test that is one o£ chi- 
square goodness of fit test, Kolmogorov-Smimof goodness of fit test; liffi 
normality and Jarque-Bera test of normality. 

83. The computer program product of Oaim 82, wherein said iKnxnaldi^^ 
one of a nonnal cumulative distribution function and a nonnal probability distribution 

« 

function in accordance with said at least one normality test 

84. The computer program product of Qaim 83, wherein said distribution of 
observed data associated with a component approximates one of: a Gaussian distribution if 
said component is healthy and a non-Gaiissian distribution otherwise. 



110 



WO02/09Sd33 PCTAJS02/16380 

85. The computer program product of Claim 82, further comprismg machine 
executable! code for: 

determining a number of differences between said observed data and eT^ected data, 
said expected data being represented by said normal distributicm; and 
detennining a sum using said differences; and 

detennining &at said observed data is not nonnally distribtited, said CTtical 
being a threshold if said number of differences is greater than a critical vahie. 

86. The computer program product of Claim 85, further comprising:* 

machine executable code for determining a score being a maximum deviation finom 
said critical value, said condition indicator being said scene. 

87. The computer program product of Claim 86, wherein sensitivity of said condition 
indicator increases as a number of observed data Vahies increases. 

» 

88. The computer program product of Qaim 81, wherein said normal distribution 
approximates a distribution of expected values. 
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89. A method executed in a computer system for detemiining a conditioii indicator 
associated witb a component, the method comprising: 

determining a total impulse signal in accordance with configuration data, said total 
impulse signal being a superposition of gear and bearing noise represented as a convolutioii 
of a gear and bearing signal Avith a geaibox transfer function; and 

determining a condition indicator in accordance with said total impulse sigoaL 

90. The method of Claim 89, further comprising: 

representing a total impulse signal generated by a configuration of associated with 
said conq>onent as: 

[impulse]® f(Gear)iS>f (Bearing)^ /{Case}^ [impulse]® [f (Gear)® /(Bearing)® /{Oise)] 

inwhidi ® represents a convolution operation. 

91. The method of Claim 90, furthar comprismg: 

representing convolution operations in a time domain to equivalent opraiations in a 
firequency domain. 

92. The method of Claim 90, further comprising: 

estimating [/{Gear) 0 /(Bearing)® /(CaseJ] as a transfer function in a fi?equeiicy 
domain using a linear predictive coding technique to deconvolute a signal into its base 
components. 



93. The method of Claim 92, further comprising: 
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estiiiiating said transfer function, H, in said frequency domain as a/B, wherein a«= 
(al, an), each ai representing an ith coefficient for an orderp, n=p+l, as: 

and B is an estimate of an error represented as: 
in which: 

b (y-yte)^ y = y[i. 2, n], 

■ . ■ . I • 

Yh0 is an estimated value of y, jkat = ax, 
X is a time delayed signal represented as: 



a^n— l,n— 2,.ji— p] 
- 2,11 - 3,n - p -1] 



where a = (x^xX' x^y, values fin-al an. 



94. The method of CHaim 93/iurther comprising: 

estimating an impulse, IMP, in said fiiequency domain of said componexit as: 

IMP = exp(log(Y) - log(H)), 
inT^ch: 

Y = ffl(y) and H = fR(h), where fft is the Fourier transfomi function, y and h are in a 
time domain, Y and H are in said frequency domain. 

95. The method of Claim 94, wherein a value associated with H increases as a fault 
increases. 
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96. The method of Claim 94, wherein said condition indicator is said value of IMP. 

97. The method ofClaim 94, further comprising: 

calculating a power spectral density of said impulse IMP in a fiequency domain; and 

deteniiininga value of said power spectral density at a frequency of interest, s ' 
condition indicator being said value. « 

98. The method of Claim 97, wherein said frequency of interest is at least one of: a 
bearing passing freqiiency for a bearing fault, and a mesh frequency for a gear fiuilL . 

99. The method of Claim 98, fijrflier comprisiiig: 

performing a Fourier transibxmation to obtain IMP in said fi-equency domain. 

100. The method of Claim 89, fbther comprising: 

detecting a fault in coimection with predetemiined values of said healHi status iigmg 
said condition indicator, wherein said faiSlt being detected is one of a pit and spaU on one o£ 
a gear tooth, ioner bearing race, outer bearing race^ and bearing roller element 
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101. A computer program product for deteirnining a conditi^^ 
with a component, the computer program product comprising madim 

determining a total impulse signal in accordiance with configuration data, said total 
5 impulse signal being a supeipositipn of gear and bearing noise represented as a conyolutioii 
of a gear and bearing signal with a geaibox transfer function; and^ 

determining a condition indicator in accordance with said total inqjulse dgntiL 

1 02. The computer program product of Claim 1 01 , further xx)niprisiii^ machiiie 
10 executable code fan 

representing a total impulse signal generated by a configuration of associated wifli 
said component as: 

\impulse]^ /(Gear)^ /(Bearing)® f {Case) s [impulse](B> [f{Cear)® flBearhtg)^ f{Case)] 

15 

in which 0 represmtsaconvoluttoa operation. 

103. The computer program of Claim 102, further comprising machine executable 
code for: 

20 representing convolution operations in a time domain to equivalent <q>erations in a 

frequency domain* 
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1 04rThe computer program product of a^m 102, furihercomp . 
executable code for 

estimating [/"(Gear)® /(Bearing)® /(Gwc)] as a transfer flmction in a fiequency 

r 

domain using a linear predictive coding technique to deconvolute a signal into its base 
S components* 

105. The computer program product of Claim IM/fiirther comprising m 
executable code fon 

estimating said transfer function, in said frequency domain as a^, herein a « 

10 (al»«., an), each d representing an ith coefficient for an order p,n»p+l,^w 

= - 1]+ a^^ii - 2]+ a3x[ii - 3]+ . . . 

and B is an estimate of an error represented as: 
15 B = 2;flffb 
in which: 

b « 0^.yA«)^ y « y£l, 2, ... n], 
20 y^at is an estimated value of y, = ax, 
X is a time delayed signal repr^ented as: 



x= 



'x[n — 1, n - 2,. ji — jo] 
jc[n - 2, n - 3, II - j> - 1] 



wherBa = (x*^x)** x^y, values for al .. an. 

25 

1 06. The computer program product of Qaim 1 05, further comprising machine 
executable code fon 



116 



wo 02/095633 PCT/US02/16380 

estimating an impulse^ IMP, in said frequency domain of said component as:. 

IMP==exp(log(Y).log(H)X . 
in which: ' ; , 

Y = ffl(y) and H = fRQx). where fft is the Fonriw transform function, y and h are in a 
time domain, Y and H are in &aid frequency doniai^ 

107. The computer program product of Claim 106, wh^ein a value associated witii H 
increases as a fault inbeases. 

1 08. The computer program product of Qaim 1 06, wherein said condition indicator 
is said value of IMP. /" 

1 09. The computer program product of Claim 106, fiirther comprisiiig madbine 
executable code for: 

calculating a power spectra] dendty of said impidse IMP in a fiequency domain;^ a 
determining a value of said power spectral density at a frequency of interest, said 
condition indicator being said value. 

110. The computer program product of Claim 1 09, wherein said frequency of intoest 
is at least one of: a bearing passing frequency for a bearing faul^ and a mesh frequency for a 
gear fault 

111. The computCT program product of Claim 110, fiuther comprising machine 
executable code for . 

performing a Fourier transformation to obtain IMP in said frequercy domain. 
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112. The computer program product of Qaim 111, further comprising 
executable code fan 

detecting a fault in connection with predetennined values of said health status using 
said condition indicator, wherein said fault being detected is one of a pit and spall on one of: ' 
a gear tooth, inner bearing race, outer bearing race, and beating roller dlemeot 
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. 113. A method executed in a computer syrtem for detennin^ 

component udng at least one condition indicator, flie method comprisiiig: 
detennining said at least one cbndition indicator usii^ at least one of: ap impulse 
determination technique and a statistical normality tes^ and 

determiningsaidhealthindicatorinaccordance with said at least one conditi 
indicator. * 

1 14* The method of Claim 1 13, wherein said statistical normality test is one of: chi- 
square goodness of fit test, Kblmogorov-Smimof goodness of fit test, lilliefcnrs test of 
noimality and Jaique-Bera test of normality. 

115. Tie method of Claim 113, wherein expected data values approximate a normal 
distribution. 
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116. A computer program product for detemiining a health status of a component 
using at least one condition indicator, the method comprising: 

determining said at least one condition indicator using at least one o£ an impulse 
determination teclmique and a statistical normality test; and 

determining said health indicator in accordance with said at least one conditian 
indicator. 

117. The computerprogramproduct of Qaim 116, wherein said statistical noimalify 
test is one of: chi-square goodness of fit test, Kolmogorov-Smimof goodness of fit test, 
Lilliefors test of normality and Jarque-Bera test of normality. 

118. The con^uter program product of Claim 116, wherein expected data values 
approximate a normal distributioa. 
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